METHOD AND SYSTEM FOR DETERMINING INSULIN DOSING 
SCHEDULES AND CARBOHYDRATE-TO-INSULIN RATIOS IN DIABETIC 

PATIENTS 
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RELATED APPLICATION 

This application claims priority to and the benefit of provisional applications, Serial No. 
60/456,271 filed March 19, 2003, Serial No. 60/532,487 filed December 26, 2003, and 
10 application filed February 12, 2004 (Express Mail No. ER946670647US). 

FIELD OF THE INVENTION 

The invention relates generally to the field o f d igital aids to assist in the treatment o f 
1 5 diabetic patients who use insulin pumps or multiple dosing insulin regimens and provides 
a method for determining insulin dosing schedules in diabetic patients. 

BACKGROUND OF THE INVENTION 

20 Diabetes Mellitus has been treated for many years by insulin injection. Three recent 
advances are changing diabetes care: The Insulin Pump, new insulin formulas for 
Multiple Dose Injection, and Inhaled Insulin. These are discussed below: 

THE INSULIN PUMP: The invention of the insulin pump revolutionized diabetes care. 

25 It is a battery-powered device about the size of a pager. It contains a cartridge of insulin 
and pumps the insulin through a flexible tube into the patient via an "infusion set", which 
is a small plastic needle or "canula" fitted with an adhesive patch. The invention of the 
pump makes it possible to adopt a typical insulin regimen as follows: Basal Insulin is 
injected slowly and continuously at a rate that can be programmed to change multiple 

30 times during the day (about 4 or 5 changes per day is common). Between the changes, 
the Basal Insulin Rate of infusion is constant. The constant periods are called "intervals". 
Additionally, boluses of insulin can be injected on command by the patient. There are 
two main types of boluses: 

35 Meal Boluses are infused just before a meal in an amount, proportional to the glycemic 
effect of the meal. This is generally proportional to the number of grams of 
carbohydrate in the meal. The proportionality constant is a personalized number called 
the Carbohydrate-to-insulin Ratio, CIR. It is used as follows: 

40 Meal Insulin Bolus = (grams of carbohydrates in the meal)/CIR ( 1 ) 

This calculation is generally performed by the patient, but there are pump models that 
can store the patient's CIR in memory and require only the grams of carbohydrate in 
the meal as the input. 



Correction Boluses are infused immediately after a Blood Glucose test has been 
performed; the amount of the correction bolus is proportional to the error in the blood 
glucose concentration from the patient's personalized Target Blood Glucose. The 
proportionality constant is a personalized number called the Correction Factor, CF. It is 
5 used as follows: 

Corrective Insulin Bolus = (Blood Glucose concentration - Target)/CF (2) 

There are two types of Corrective Bolus, each with a different Target: 

10 

Time-Boundary Corrective Insulin Boluses are administered in a fasting state at the end 
of a time interval. 

After-Meal C orrective Boluses are administered from one to five hours after a meal, 
15 most often within the time interval. 

Recently, p ump m anufacturers h ave b een i ncorporating d igital features in their pumps 
that make treatment easier. Some pumps can store the values of CF and Target and 
require only the Blood Glucose Concentration (BG) as input. Among these new digital 
20 features is the "Insulin-On-Board" feature. This feature mathematically models the 
amount of insulin still in the body at a given time after a bolus and recommends 
reductions to the boluses accordingly. This feature makes After-Meal Corrective Boluses 
more safe and practical. 

25 MULTIPLE DOSE INJECTION (MDI): Advances are being made in developing 
different types of insulin. Some are very long acting and non-peaking. The long-acting 
insulin can be injected as infrequently a s o nee p er d ay i n a regimen v ery s imilar t o a 
pump patient's basal insulin regimen. Injections of rapid-acting types of insulin can be 
given as meal and correction boluses. The two types together act as a system. These 

30 insulins are available in portable "pens" (named for their resemblance to writing 
implements). The pens have been mated with BG meters in "kits" in which the devices 
communicate so that the combined memory is stored in one of the two devices in the 
"kit". 

35 INHALED INSULIN: Inhaled insulin delivery systems are under development for short- 
acting i nsulin. It i s e xpected t hat t he i nhalers w ill b e c ombined with BG meters into 
"kits" like the ones used for MDI, then the present invention will be able to handle 
inhaled insulin in the same manner. This development is expected in the future. 

40 The nature of diabetes care is very quantitative. Ironically, the proliferation of numbers 
makes the use of lengthy algorithms on pocket calculators too time-consuming and 
therefore prohibitively expensive. The majority of endocrinologists, therefore, use 
experience-based subjective methods. In the interest of providing greater subjective feel 
for the case at hand, endocrinologists often use the numbers to simply help them discern 
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trends; then treat the trends. For instance, they commonly view the blood glucose (BG) 
scatter charts and printouts to discern trends by such subjective means as the visual 
density of dots on the BG scatter chart and the relative location of the areas of highest 
density. They translate these trends into insulin dose changes using their experience. 

5 

Experienced-based and subjective methods are often not uniform from one practitioner to 
the next. Additionally, there is a shortage of endocrinologists and other diabetes 
specialists. Accordingly, the management of diabetes is done in a disorganized manner by 
clinicians of widely varying degrees of expertise. The result is that control of diabetes in 
10 most patients, while satisfactory, is not optimal. As a result of sub-optimal BG control, 
the course of diabetes can include complications involving all body systems. These 
complications are associated with premature mortality and are associated with a cost, 
which amounts to 19% of the health dollars to care for 6% of the population. 

15 BRIEF SUMMARY OF THE INVENTION 

The present invention overcomes the aforementioned disadvantages of current care, by 
providing a method to analyze and prescribe changes to the daily insulin-dosing schedule 
of diabetic patients using insulin pumps, multiple-dose subcutaneous injection, or inhaled 

20 insulin. The method divides the patient's day into selected time intervals in which 
adjustable schedules are provided for Basal Insulin dosage rates and Carbohydrate-to- 
Insulin Ratio(s) (to determine Meal Insulin doses). The time boundaries and Basal Rate 
changes are usually set by the Practitioner to coincide with the patient's meals. The 
patients are usually encouraged to test their Blood Glucose (BG) at the time boundaries 

25 just before they eat. The invention uses two systems of insulin nomenclature: 

CONVENTIONAL INSULIN NOMENCLATURE: 

Total Daily Insulin = Prescription Insulin + Corrective Insulin (3) 
30 where: Prescription Insulin = Basal Insulin + Meal Insulin (4) 

The invention incorporates the concept of utilizing the Corrective Insulin over a selected 
time interval as an "error" in the patient's Prescription Insulin (Basal Insulin + Meal 
Insulin) for the whole day as well as for each time interval. Methods are included for 
35 estimating the change to one of the two components of Prescription Insulin, then 
determining the change to the other component by subtracting from the error. Therefore, 
there are two basic algorithmic forms, which are called "Floats": 

The Meallns Float: 

40 

The change in Basal Insulin is estimated first; then the invention calculates the change in 
Meal Insulin to be the error in Prescription Insulin minus the change in Basal Insulin. 
Some of the ways for estimating the change in Basal Insulin are: 
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Borrowing Basal Rate from another interval; 
Estimating Basal Rate from another algorithm (e.g. a Basal Float); 
Estimating change in Basal Insulin from the carb-free latter part of a time interval 
using the Time-Boundary Corrective Insulin (at the end of the interval) as the 
5 error indicator. The float is done on the first part of the interval, using the After- 

Meal Corective Insulin as the error 

or 

The Basal Float: 

10 The change in Meal Insulin is estimated first; then the invention calculates the change in 
Basal Insulin to be the error in Prescription Insulin minus the change in Meal Insulin. 
Some of the ways for estimating the change in Meal Insulin are: 

Borrowing Meal Insulin or CIR from another interval; 
15 Estimating the change in Meal Insulin for an interval as a share of the Change in 

Total Day's Meal Insulin in the same proportion as (Carbs for the interval)/(Total 
Day's Carbs); 

Estimating the change in Meal Insulin for an interval as a share of the Change in 
Total Day's Meal Insulin in the same proportion as (Meal Insulin for the 
20 interval)/(Total Day's Meal Insulin). 

ENHANCED INSULIN NOMENCLATURE: 

Also, to provide another way of accounting for After-Meal Insulin, the invention 
identifies "Enhanced Variables" as a system of insulin nomenclature that lumps the 
25 After-Meal Corrective Insulin with the Meal Insulin as follows: 

Total Daily Insulin=Enhanced Prescription Insulin+Time-Boundary Corrective Insulin (5) 

where 

30 

Enhanced Prescription Insulin = Basal Insulin + Enhanced Meal Insulin (6) 
and: 

35 Enhanced Meal Insulin = Meal Insulin + After-Meal Corrective Insulin (7) 
The Floats are very similar: 
The Meallns Float: 

40 

The change in Basal Insulin is estimated first; then the invention calculates the change in 
Enhanced Meal Insulin to be the error in Enhanced Prescription Insulin minus the change 
in Basal Insulin. Some of the ways for estimating Basal Insulin are: 
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Borrowing Basal Rate from another interval 

Estimating Basal Rate from another algorithm (e.g. a Basal Float) 



or 

5 

The Basal Float: 

The change in Enhanced Meal Insulin is estimated first; then the invention calculates the 
change in Basal Insulin to be the error in Enhanced Prescription Insulin minus the change 
10 in Enhanced Meal Insulin. Some of the ways for estimating Enhanced Meal Insulin are: 

Borrowing Enhanced Meal Insulin or CIR from another interval; 

Estimating the change in Enhanced Meal Insulin for an interval as a share of the 
Change in Total Day's Enhanced Meal Insulin in the same proportion as (Carbs 
1 5 for the interval)/(Total Day' s Carbs); 

Estimating the change in Enhanced Meal Insulin for an interval as a share of the 
Change Total Day's Enhanced Meal Insulin in the same proportion as 
(Enhanced Meal Insulin for the interval)/(Total Day's Enhanced Meal Insulin); 

20 

IN GENERAL: 

It may not be desirable to apply all of the error term to effect a change. Methods are 
included for applying a limited amount of the error in order to avoid overshoots. 

25 The above selection of "Float" algorithms have been developed in two main versions: 

Daily Update Version: The Float algorithms above are embodied in a program that uses 
the previous day's data to make calculations for the present day. In this version, the full 
amounts of the error terms are not applied, but instead a predetermined fraction of the 
30 error terms are applied. If the error was valid, a reduced amount of it will show-up the 
next day. This next day's reduced error will be reduced by the same fraction and so on 
until the error disappears asymptotically. This version is especially well-suited for 
installation in insulin pumps, inhaled insulin kits, MDI injection kits, PDA's and other 
portable devices. 

35 

Multiple Days' Data Version: The Float algorithms above are embodied in a program 
that uses the recent calendar period (e.g. last few weeks) for the data source. The data are 
therefore averages. The full amounts of the error terms are not necessarily applied, but 
instead the fraction or amount of error to be applied is input by the practitioner. These 
40 versions are suitable for installation in the practitioner's computer. A sub- version of this 
version determines the error fractions automatically. It is suitable for the patients' 
computer, or a website, for use with patients who have insulin pumps, or kits for Multiple 
Daily Injection or Inhaled Insulin. 
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The various aspects of the present disclosure may be more clearly understood and 
appreciated from a review of the following detailed descriptionof the disclosed 
embodiments and by reference to the appended figures. 

5 BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 illustrates an exemplary embodiment of the Input Form for the Multiple Days' 
Data Version for Insulin Pumps (Type D). The inner panel is the SubForm. It has two 
pages, which can be reached by scrolling. The first of these is shown. It contains single- 
10 valued data (not scheduled data). 

Figure 2 illustrates an exemplary embodiment of the Input Form for the Multiple Days' 
Data Version for Insulin Pumps (Type D). The inner panel is the SubForm. It has two 
pages, which can be reached by scrolling. The second of these is shown. It contains daily 
15 schedule data. 

Figure 3 an exemplary embodiment of the source of the multiplier FinsAuto, which 
automates dRxInsAuto in the automated Multiple Days' Data Version (Automatic Digital 
Advisor) by comparing the standard deviation of the patient's BGs to the mean standard 
20 deviation from the database. If the patient's number is high, the multiplier (which is 
always < 1) is reduced. The automated change in Enhanced Prescription Insulin is 
FinsAuto times the Time-Boundary Corrective Insulin times another factor. 
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DETAILED DESCRIPTION OF THE INVENTION 



1. TABLE OF CONTENTS 

The invention is a set of algorithms used to determine the various insulin-dosing rates, 
and Carbohydrate-to-insulin Ratio(s) that comprise the daily insulin-dosing schedule for 
a patient. The "Float Algorithms" apply to all exemplary embodiments described herein. 
30 The description, as shown in the table of contents below, is broken down first by 

"Embodiment" (pump or MDI) and then by "Version" (Daily Update or Multiple Day's 
Data) and finally by "Algorithm" ( Basal Float 1, Meal Insulin Float 1, etc). 



1. TABLE OF CONTENTS 

2. SUFFIXES and TIME INDICES 

3. GLOSSARY 

40 4. GENERIC DERIVATIONS 

4.1 MAKING CHANGES IN THE CARBOHDRATE-TO-INSULIN RATIO, CIR 
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4.2 SOME IMPORTANT EQUATIONS 

4.2.1 CONVENTIONAL INSULIN NOMENCLATURE 
5 4.2.2 ENHANCED INSULIN NOMENCLATURE: 
4.2.3 GENERAL 

4.3 THE TWO-LEVEL BASAL SYSTEM, and the Kf CALCULATOR 

10 

4.3.1 STATISTICAL CORRELATION 

4.3.2 AVERAGING 

15 5. TYPES OF INSULIN DELIVERY SYSTEMS 

6. DESCRIPTION BY. VERSION. EMBODIMENT. ALGORITHM 
6.1 VERSIONS USING DAILY UPDATE 

20 

6.1.1 PRELIMINARY DERIVATIONS 

6.1.1.1 TIME INTERVALS AND CORRECTIVE INSULIN 
25 6.1.1.2 Krxlnsl: GOVERNS THE SIZE OF INSULIN CHANGES 

6.1.2 FOR PUMPS (TYPE E) 

6.1.2.1 BASAL-TO-TOTAL RATIO 

30 

6.1.2.2 BASAL FLOAT 1. ( deMeallns is proportional to eMeallns or Carbs) 

6.1.2.3 BASAL FLOAT 2. ( deMeallns value from outside the interval) 
35 6.1.2.4 MEAL INSULIN FLOAT 1 

6.1.2.5 OVERVIEW OF 6.1.1.2 and 6.1.1.3 

6.1.2.6 MEAL FLOAT 2 (uses AMCorlns(i) as error ) 

40 

6.1.3 FOR MULTIPLE DOSE INJECTION (MDI) AND INHALED INSULIN: 
6.1.3.1 BASAL FLOAT 1. (not used with MDI or Inhaled Insulin') 
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6.1.3.2 BASAL FLOAT 2. ( The value of dMeallns is from outside the interval ) 

6.1.3.3 MEAL INSULIN FLOAT 1 
5 6.1.3.4 OVERVIEW OF 6.1.3 

6.1.3.5 MEAL INSULIN FLOAT 2. (uses AMCorlns(i) as an error term) 

6.1.4 SKIPPED BG'S 

10 

6.1.5 A MODIFIED BOLUS CALCULATOR 

6.1.6 CHANGING THE PATIENT'S SCHEDULE 
15 6.2 VERSIONS USING MULTIPLE-DAYS ' DATA 

6.2.1 MULTIPLE DAYS' DATA (DIGITAL ADVISOR) FOR PRACTITIONERS 
6.2.1.1 FOR PUMPS 

20 

6.2.1.1.1 PUMP TYPE D and E 

6.2.1.1.1.1 Basal Float 1 , (deMeallns is proportional to eMeallns or Carbs) 
25 6.2. 1 . 1 . 1 .2 Basal Float 2, (not used with Multiple Days' Data) 

6.2.1.1.1.3 Meal Insulin Float 1 

6.2.1.1.1.4 Overview of 6.2. 1 . 1 . 1 . 1 and 6.2. 1 . 1 . 1 .3 Multiple Days' Data Type D 
30 pumps 

6.2.1.1.1.5 Meal Insulin Float 2 (uses AMCorlns as the error ) 

6.2.1.1.2 PUMP TYPE C 

35 

6.2.1.1.3 PUMP TYPE B 

6.2.1.1.4 PUMP TYPE A 

40 6.2.1.2 SUBCUTANEOUS MULTIPLE DOSE and INHALED INSULIN 
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6.2.1.2.1 BASAL FLOAT 1 

6.2.1.2.2 MEAL INSULIN FLOAT 1 

5 6.2. 1.2.3 MEAL INSULIN FLOAT 2 (uses AMCorlns(i) as an error ) 

6.2.1.3 LIMITED DOMAIN OF deRxInsl and the SAFETY-NET FORMULA for 
MULTIPLE DAYS' DATA VERSIONS 

10 6.2. 1 .3. 1 LIMITED DOMAIN OF deRxInsl 

6.2.1.3.2 SAFETY-NET FORMULA 

6.2.1.3.2.1 For Basal Float 1 

15 

6.2.1.3.2.2 For Meal Insulin Float 1 

6.2.1.3.2.3 For Meal Insulin Float 2 (uses AMCorlns(i) as error ) 
20 6.2.2 AUTOMATIC MULTIPLE DAYS' DATA (DIGITAL ADVISOR) 

6.2.2.1 AUTOMATION OF deRxInsl 

6.2.2.2 AUTOMATION OF deMeallns 

25 

6.2.2.3 AUTOMATIC ROUNDING OF CIR FOR PUMP TYPE A 
6.3 NON-FLOAT ALGORITHM 

30 7 TIME INTERVALS 

7.1 FINDING THE DAY'S LAST ENTRY 

35 2. SUFFIXES and TIME INDICES: 

In the present invention, the 24 hour day is divisible into multiple time intervals. The 
boundaries of time intervals are numbered in the manner of: timeO, timel, time2, etc. 
These time boundaries will be referred to generically as time(i), where "i" denotes "time 
40 index". 

The Time Intervals, dt(i), are found by subtracting the consecutive time boundaries. 
dt(i) = time(i+l)-time(i) (8) 
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Each time interval is numbered the same as its upstream time boundary as shown above. 
Special mention should be made of the interval surrounding midnight. 



5 dtO = 24 + timel - Tmax (9) 

where Tmax is the last time boundary in the day, usually Bedtime. The method of 
finding Tmax is set forth in a later section. 

10 The time boundaries are generally defined by the changes in the Basal Rate (in pumps) or 
by average mealtimes. Practitioners often set up the patient's schedule so that the Basal 
Rate changes, mealtimes, BG tests, Corrective Insulin Boluses, and Meal Insulin Boluses 
all occur at these time boundaries. While not mandatory, this practice is preferred for use 
with the invention. The Basal Rates are numbered the same as time intervals, thus BR1 

15 goes between timel and time2. 

The other parameters (besides Basal Rate) are identified with the exact time boundaries, 
though in practical use, the patients cannot be expected to follow such a strict schedule. 
Therefore, there is a need for another set of time periods (called "bins" to distinguish 

20 them from intervals), each of which envelops a time boundary. The parameters that are 
said to occur "at a time boundary" actually occur in a bin and are given the same index 
number as the time boundary enveloped by the bin. The tests and boluses can be sorted 
into these "Time-Boundary bins". In the case of the Multiple Days' Data Versions, this 
is most often handled by the downloading software of the Blood Glucose Meters and 

25 Insulin Pumps, in which a schedule of bins can be set up. In the case of the Daily Update 
Versions, the midpoints of the primary time intervals may be used to define the 
boundaries of these bins. 

Note that Time-Boundary Corrective Insulin measures the insulin error in the previous 
30 interval, thus when calculating parameters for the ith interval, the insulin error is 
"TBCorIns(d,i+l)". 

Parameters that do not follow the schedules of time intervals lack the time interval index 
integer. In the Multiple Days' Versions, these parameters lack any index at all. 

35 

In the Daily Update Versions, all parameters are given a day's index, "d", "d-1", etc in 
the manner of BF(d,i) or TDD(d). 

Figure 1 and Figure 2 show the two pages of an exemplary Input Form for Version 
40 6.2.1.1.1, the Multiple Days' Data Version for Pump D. Figure 1 shows un-indexed 
parameters that have a single value for the Interaction. Figure 2 shows the time intervals 
containing time-interval-indexed parameters. 
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The parameters that pertain to a unique date and time are denoted with the functional 
form: "parameter(t)". A few of these appear in the derivations herein for explanatory 
purposes, but none are used in the present invention. 

5 3. GLOSSARY: 

Average Glycemic Index (AGI): For average carbohydrates, the ratio of the grams of 
glucose entering the blood within two hours divided by the total grams of sugar in the 
carbohydrate. There are published studies for this figure. It ranges between 60% and 
10 90%. 
s 

AIM Formulas : An acronym for "Accurate Insulin Management" (See Ref 1, which is 
incorporated by reference as if fully set forth herein). A set of statistically derived 
formulas used to estimate three parameters for the patient. There are three formulas. In 
15 each, a constant (beginning with "K") was determined statistically: 

BasalAIM = Kb * TDDavg (10) 
CFaim = Kef / TDDavg (1 1) 

CIRaim = Kcir * BodyWeight(lbs) / TDDavg (12) 

20 

The values of these constants change from time to time as new research is done. 
The latest values (See Ref 1) are: Kb = 0.48, KcfH700, Kcir=2.8 

After-Meal Corrective Insulin: Corrective boluses taken after meals (post-prandially) in 
25 the previous day in the ith interval. (Daily Update Version) 

AMCorIns(d-l,i): "After-Meal Corrective Insulin ": Corrective boluses taken after meals 
(post-prandially) in the previous day in the ith interval. (Daily Update Version) 

30 AMCorlns(i): "After-Meal Corrective Insulin for the ith time interval" The average of 
corrective boluses taken after meals (post-prandially) in the ith time interval. (Multiple 
Past Days Version) 

Basal Rate: The rate in (insulin units / hour) at which Basal Insulin is administered. 

35 

Basal Rate Float: CIR or Meal Insulin is determined by the practitioner or by an 
estimation formula and BR is determined by subtracting the change in Meal Insulin 
(Enhanced or conventional) from the total desired change, dRxInsl (Enhanced or 
conventional). 

40 

Basal(d): "Basal Insulin." Basal insulin for the current day, calculated from the basal 
rates by the invention. (Daily Update Version) 

Basal(d,i): Basal Insulin for the current day and time interval (Daily Update Version) 
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Basal(i): Basal for the time interval (Multiple Days' Version) 

Basal: "Basal Insulin." 1. General definition: Insulin that is administered continuously 
5 by pump or insulin that is injected manually and stays in the body for a long time due to 
its special chemical make-up. 2. Computer variable: The current total of Basal I nsulin 
administered during the day as calculated from the basal rates. 
SUM( Basal(i)) over the day calculated by the invention. (Multiple Days' Version) 

10 BasalTot: "Basal Insulin Total." Current Basal Insulin dose programmed into the pump 
at the time of the interaction with the Practitioner. Calculated by the pump. It should 
equal Basal. (Multiple Day's Version) 

BG(i): Average of blood glucose tests over many days at approximately the same time 
15 boundary, time(i), during the day, calculated by the glucose meter's downloading 
software. (Multiple Day's Version) 

BGmean: Overall average of blood glucose tests since last interaction or in the calendar 
period being analyzed, obtained from the BG meter's downloading software. (Multiple 
20 Day's Version) 

BGpd: "BG's per day." The average number of BG tests per day in the calendar period 
being analyzed. (Multiple Day's Version) 

25 BGsd(i): "BG standard deviation for the ith time interval" obtained from the BG meter's 
downloading software. (Multiple Day's Version) 

BGsd: "BG standard deviation," over the calendar period being studied. Obtained from 
the BG meter's downloading software. (Multiple Day's Version) 

30 

BG, Blood Glucose Concentration: These tests are self-administered several times during 
the day, usually on the time boundaries. At the present state of the art, they are 
determined by placing a drop of blood on a test strip in a blood glucose meter. The meter 
stores the values and times in its memory. At each practitioner/patient interaction, the 
35 values are generally downloaded from the Blood Glucose Meter's memory by software 
that graphs the data points and calculates s everal values including the overall average 
denoted BGmean, and the time boundary averages, denoted BG(i). Some models of BG 
meter are able to link to insulin pumps providing a combined and integrated download of 
data. 

40 

Bin: A time period enveloping a Time Boundary of the primary Time Interval system. 
Any BG test or bolus occurring within the Bin is treated as though it occurred on the 
Time Boundary. 
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Bolus: {from the Latin, ball) Insulin infused in a short elapsed time as ordered by the 
patient, as distinguished from Basal Insulin, which is infused slowly, continuously and 
automatically either by a pre-programmed pump or by injection of a slow-acting insulin. 

5 Bolus Time Period(i): Boluses are identified with Time Boundaries, but there is 
variability in the timing of them. To resolve this, a system of Bolus Time Periods is used. 
It has boundaries alternating in staggered fashion between the boundaries of the regular 
Basal Intervals such that any BG or Bolus falling in the Bolus Time period is 
automatically considered as occurring exactly on the regular Basal Time Boundary. 

10 

Boli(i): "Bolus Insulin": The average sum of Meal Insulin + Corrective Insulin for a 
time interval, so called because it is administered in boluses rather than continuously. 
(Multiple Past Days' Version) 

15 BoT : "Basal over Total", Basal / TDD . 

BoT(d): "Basal over Total" Basal/TDD for the day. 

BoTFbk: "Basal over Total, Feedback". A parameter in the Daily Update Versions. 
20 This factor, multiplied times the day's proposed insulin change, will yield the amount of 
change to Basal. 

BoTTgt: "Basal over Total, Target " The desired value of Basal/TDD (Daily Update 
Version) 

25 

BoTTgtRec: "Basal over Total, Target Recommended" Recommended value of BoTTgt 
(Daily Update Version) 

BR(d,i): "Basal Rate" for the ith time interval." Current setting for the rate at which 
30 Basal Insulin is delivered. The rate is programmable in the pump so that it can change 
several times during the day. (Daily Update Version) 

BR(i): "Basal Rate" for the time interval." Current setting for the rate at which Basal 
Insulin is delivered. The rate is programmable in the pump so that it can change several 
35 times during the day. (Multiple Days' Version) 

BRavg: "Basal Rate Average", Basal/24 

BRavgRec: "Basal Rate Average Recommended" 

40 

BRf: "Basal Rate, Fasting". The BR that would sustain a patient's blood glucose in the 
target range if no meals are consumed within 12 hours previous and if the measurement is 
not being made in the Late-Sleep Interval (to avoid the complications of the Pre-Dawn 
effect). 
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BRkey: The key Basal Rate. Other Basal Rates are pegged to it. 
BRlateSlp: "Basal Rate in Late-Sleep time interval". 

5 

BRrec(i): "Basal Rate, Recommended for the ith time interval." Calculated by the 
present invention. (Multiple Past Days Version) 

BRreliable: A set of basal rate data chosen for reliability, high frequency of associated 
10 BG tests, etc. Used with the Kf calculator 

BRrx(i): "Basal Rate, Prescription for the ith time interval." The value the Practitioner 
gives to the patient, after reviewing the case. (Multiple Past Days Version) 

15 BRsimilar: A set of basal rate data chosen for similarity or closeness to the fasting basal 
rate. Used with the Kf calculator 

Carbohydrate-to-Insulin Ratio ( see CIR): A personalized conversion constant; current 
value. 

20 

Carb-Counting: The technique of determining meal insulin by first counting the grams 
of carbohydrate in the meal about to be consumed and then dividing by CIR. 

Carbm: "Carbohydrate grams modified"; Used in bolus calculators that incorporate 
25 exercise. 

Carbs: "Grams of Carbohydrate" 

CarbSh(i): "Carb Share." When used with all but pump type A, actual grams of 
30 carbohydrate or 15 gm "exchanges" are used. When used with pump type A, with no 
memory for Meallns(i), this parameter is the Practitioner's or Patients's estimate of the 
relative magnitude of carbohydrates consumed in each time interval. For pump type A 
any units are acceptable, because t he p arameter i s o nly u sed as a p ercent o f t he t otal, 
CarbShTot. Example units: grams of carbohydrate (preferred), percent of the total, units 
35 of meal insulin (if CIR constant). (Multiple Days' Version) 

CarbShTot: "Carb ShareTotal"; calculated as CarbShTot = sum(CarbSh(i)) 

CF: "Correction Factor." A personalized factor used to calculate Correction Boluses. 
40 See CorBol(t) and AIM 

Change to: as in "Change to Pescription Insulin", "Change to Basal Insulin", "Change to 
Enhanced Meal Insulin". This phrase denotes a proposed change in the named quantity, 
from the last measured value to the recommended value. For the Daily Update version, 
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this is the change from one day to the next. For theMultiple Days 5 Data Version this is 
the change from the average over the past calendar period to the recommended value. 
Digital variables denoting changes are preceeded by a lower-case "d" as in calculus. 

5 CIR(d,i): "Carbohydrate-to-insulin Ratio" The CIR for a time interval dt(i) in the current 
day. (Daily Update Version) 

CIR(i): "Carbohydrate-to-insulin Ratio" The CIR for a time interval dt(i). (Multiple Past 
Days Version) 

10 

CIR: "Carbohydrate-to-insulin Ratio". A personalized conversion constant; current 
value. 

CIR = (weight of carbohydrates consumed) / (insulin required to metabolize the 
1 5 carbohydrates) ( 1 3) 

See MealBol(t), and AIM. Calculation of this parameter is one of the chief goals of 
PumpMaster. (Multiple Past Days Version) 

20 CIRrec(i): f, Carbohydrate-to-Insulin Ratio, recommended for a time interval" (Multiple 
Past Days' Version) 

CIRrec:"Carbohydrate-to-Insulin Ratio, recommended" (Multiple Past Days Version) 

25 CorBol(t): "Corrective Bolus at a unique date and time." The patient periodically tests 
his/her Blood Glucose concentration. If it is too high, the patient calculates a bolus of 
insulin by the formula: 

CorBol(t) = (BG(t) - Target) / CF (14) 

30 

Where Target is the desired BG level and CF is a personalized Correction Factor. 

Corlns(i): "Corrective Insulin." The sum of CorBol(t) at all times within the interval, 
including "After-Meal" and "Time-Boundary" times. (Multiple Past Days' Version) 

35 

CorIns(d,i) "Corrective Insulin." The sum of CorBol(t) at all times within the interval 
including "After-Meal" and "Time-Boundary" times during the current day. (Daily 
Update Version) 

40 CorTot(d): "Corrective Insulin Total." The sum of corrective insulin at all times within 
the current day. (Daily Update Version) 

CorTot: "Corrective Insulin Total." The average of the sum of corrective insulin at all 
times within day. (Multiple Past Day's Version) 
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dBasalAuto: The automated version of dBasal. (Automatic Multiple Days' Version) 
dBaslToAIM: The change in Basal necessary to achieve the BasalAIM value 

5 

deMealIns(d,i): "Change in Enhanced Meal Insulin" for the time interval for the current 
day. (Daily Update Version) 

deMeallns(i): "Change in Enhanced Meal Insulin" for the time interval. (Multiple Past 
10 Days Version) 

deMeallns: "Change in Total Enhanced Meal Insulin." For the whole day. (Multiple 
Past Days' Version) 

15 deMeallns Auto: The automated version of dEmeallns. (Automatic Multiple Days' 
Version) 

deRxInsAuto: "Change in Enhanced Prescription Insulin, Automatic"; estimated 
automatically. (Automatic Multiple Days' Version) 

20 

deRxInsl(d,i): "change in Enhanced Prescription Insulin in a time interval of the current 
day" (Daily Update Version) 

deRxInsl(i): "change in Enhanced Prescription Insulin in a time interval" (Multiple Past 
25 Days Version) 

deRxInsl: " change i n T otal Day's Enhanced Prescription Insulin" (Multiple Past Days 
Version) 

30 deRxInslld(i): "change in Enhanced Prescription Insulin, large domain in the interval" 
The "Safety-Net" formula. One of two formulas for deRxInsl(i). (Multiple Days' 
Version) 

deRxInslsd(i): "change in Total Enhanced Prescription Insulin, small domain"; One of 
35 two formulas for 

deRxInsl(i) for "Safety-Net" purposes. (Multiple Days' Version) 

dt(i): "difference in time"; The length of the "ith" time interval 

40 

Early-Sleep Time Interval: The time interval starting with Bedtime and ending with the 
Mid-Sleep time boundary. Most diabetic patients are encouraged not to snack at bedtime 
and to take a mid-sleep BG at about 3:00 AM. If this advice is followed, then the Early 
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Sleep Time Interval gives a good indication of the true "inactive basal rate", BRf, that 
underlies many other time intervals. 

eMealIns(d,i): "Enhanced Meal Insulin for the time interval." The sum of the Meal 
5 Boluses plus the After-Meal Corrective Boluses during time interval dt(i) for the current 
day. (Daily Update Version) 

eMeallns(i): "Enhanced Meal Insulin for the time interval." A multiple-day average of 
the sum of After-Meal Insulin plus Time-Boundary Corrective Insulin in a time interval. 
10 (Multiple Past Days Version) 

eMealInsRec(d,i): "Enhanced Meal Insulin Recommended for the current day and time 
interval". (Daily Update Version) 

15 eMeallnsRec(i): "Enhanced Meal Insulin for the time interval, recommended". 
Calculated by the invention. (Multiple Past Days Version) 

eMeallnsTot(d): "Enhanced Meal Insulin Total for the current day" The sum of 
eMeallns(i) during the current day. (Daily Update Version) 

20 

eMeallnsTot: "Enhanced Meal Insulin Total." The sum of Meallns(i) plus the sum of 
After-Meal Corrective Insulin during the whole day averaged over several days. 
(Multiple Days' Version) 

25 Enhanced: A word used to describe a certain modified variable system, in which After- 
Meal Corrective Insulin is incorporated as part of Meal Insulin. 

eRxInsl: "Total Day's Enhanced Prescription Insulin." A subset of TDD, defined as: 
The sum of Enhanced Meal Insulin + Basal Insulin during the day. It is usually 
30 "Prescription" by the practitioner. The invention proposes methods of calculating it 
automatically. 

ExerCarbs: "Exercise Carbs." The equivalent of exercise in grams of carbohydrate. 

35 FInsAuto: A multiplier that limits dRxInsAuto. (Automatic Multiple Days' Version) 

Fixed-Meal: Another simpler technique of determining meal insulin other than carb- 
counting. It involves eating meals of fixed menu preceded by fixed amounts of insulin. 

40 Float: An algorithm in which the change to one of the two main parameters, Emeallns or 
BR is estimated or determined by the practitioner, and the change to the other parameter 
is determined by subtracting the estimated insulin change from the total desired change in 
insulin. The "floated" parameter is the one that is determined by subtraction. 
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Glycemiclndex: The ratio of calories available from carbohydrates within 2 hours after 
consumption to total calories in the carbs. It is different for different types of carbs. 



Interaction: "Patient/practitioner Interaction." An inclusive term, including person-to- 
5 person interactions, and telecommunication-based interactions in which the patient's 
parameters are re-adjusted. 

Interval: The day is divided into time intervals for the purpose of enabling a patient's 
insulin dosing to be varied throughout the day and analyzed on a time-dependent basis. 
10 The intervals usually are bounded at the times of the Basal Rate changes. 

Kb: A statistically-derived constant used in the formula for estimated Basal. The latest 
publication of the "AIM" formulas give the latest value of Kb=0.48 (Ref 1). It is used in 
the formula: 

15 

Basal = Kb * TDDavg (15) 

Kef: A statistically-derived constant used in the formula for estimated CF. The latest 
publication of the "AIM" formulas (Ref 1) give the latest value of Kcir=1700 . It is used 
20 in the formula: 

CF = Kef / TDDavg (16) 

Kcir: A statistically-derived constant used in the formula for estimated CIR. The latest 
25 publication of the "AIM" formulas (Ref 1) give the latest value of Kcir=2.8 . It is used in 
the formula: 

CIR = Kcir*BodyWeight / TDDavg (17) 

30 KcirW: A constant used in the formula for estimated CER. The latest publication (Ref 1) 
give the latest value of Kcir=500 . It is used in the formula: 

CIR = KcirW / TDDavg (1 8) 

35 Kcycle: A constant representing the fraction by which CorTot is intended to be reduced 
in a specified number of days (Ncycle). (Daily Update Version) 

Keyb: as in "Key Basal Rate", An identifier for an interval that is considered to have 
dependable data and whose BR is also used in other intervals. 

40 

Keyc: as in "Key CIR", An identifier for an interval that is considered to have 
dependable data and whose CER is also used in other intervals. 
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Kf = "K fasting"; A constant <1. The ratio of fasting basal rate to a chosen "reliable" 
basal rate 

Kf=BRf/BRreliable (19) 

5 Kmauto: "K meals,auto"; a positive constant < =1 whose purpose is determining 
deMeallns. (Automatic Multiple Days' Version) 

Kfbk: "K Feedback", A constant in the formula for BoTFbk that adjusts the speed of 
convergence (days) of the Basal-to-Total Ratio to a target value. (Daily Update Version) 

10 

Krxlnsl: a positive constant < =1 whose purpose is limiting dRxInsl. (Daily Update 
Version) 

KrxInslMax: The maximum allowable value of deRxInsl/TBCorTot (Multiple Days' 
15 Version) 

Late-Sleep Time Interval: (same as Pre-Dawn Time Inteval) The time interval starting 
with the Mid Sleep time boundary and ending with the Breakfast time boundary. This 
time interval is characterized by the "Pre-Dawn Effect", which is a need for somewhat 
20 more insulin that would be normally expected. It is believed that this effect is linked to a 
release of growth hormone. This time interval is the most dependably carb-free of all the 
intervals. It would be the best candidate for the "inactive fasting basal rate" were it not 
for the Pre-Dawn Effect. 

25 Meal Insulin: Insulin taken in a "bolus" concurrently or just before consumption of 
carbohydrates. 

Meal Insulin Float: BR is determined by the practitioner or by an estimation formula, 
and Meal Insulin (enhanced or conventional) is determined by subtracting the change in 
30 Basal from the total desired change, dRxInsl (Enhanced or conventional). Then change 
in CIR is determined from the change in Meal Insulin (enhanced or conventional). Use 
of this type of Float allows the Basal Rates to be kept more uniform. This is a benefit if a 
patient skips a meal. 

35 MealBol(t): "Meal Bolus at a unique date and time." When eating, a patient calculates 
the grams of carbohydrates being consumed and calculates an insulin bolus by the 
following formula: 

MealBol(t) = ( gm of carbohydrate ) / CIR (20) 

40 

where CIR is the personalized Carbohydrate-to-Insulin Ratio 
Meallns(d,i) = Meal insulin within a time interval (Daily Update Version) 
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Meallns(i) = Meal insulin within a time interval (Multiple Past Day's Version) 

MeallnsTot = The total of meal insulin in a day 

5 NBGs(i): "number of BG's in the time interval" over the calendar period being analyzed. 
(Multiple Day's Version) 

NBGs: "number of BG's." Total number of BG tests in the BG meter since last 
download or in the calendar period being analyzed. (Multiple Day's Version) 

10 

Ncycle: the number of days until a fraction Kcycle of the insulin error is removed. 
(Daily Update Version) 

NDbg: "number of days of blood glucose." Number of days of BG tests in the calendar 
15 period being analyzed. (Multiple Day's Version) 

Peg: (verb) To set a parameter (CIR or BR) in one interval equal to a constant times the 
same parameter in another interval. The parameter is said to be "pegged" to the key 
interval. Normally the constant is equal to the ratio of the parameters on day of 
20 adjustment by the Practitioner, so that the same ratio is maintained in the future. 

Pen: A pocket-portable insulin delivery device, named for its resemblance to a writing 
implement. 

25 PmPctBGsd: The database "population" mean of the quantity (BGsd/BGmean), used in 
the automation of 

deRxInsl. (Automatic Multiple Days' Version) 

Practitioner: The physician or nurse who analyzes a diabetic patient's parameters and 
30 prescribes insulin dose regimen. 

Pre-Dawn Time Interval: (same as Late-Sleep Time Interval) 

PsdPctBGsd: The database "population" standard deviation of the quantity 
35 (BGsd/BGmean), used in the automation of deRxInsl. (Automatic Multiple Days' 
Version) 

Reliable data set: A special use of the word to indicate an interval or set of data that has 
good statistics on a certain parameter, i.e. high frequency of BG testing, low Standard 
40 Deviation. 

Rxlnsl: "Total Day's Prescription Insulin." A subset of TDD, defined as: The sum of 
Meal Insulin + Basal Insulin during the day. It is prescribed by the practitioner, unlike 
Corrective Insulin. 
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Similar data set: A special use of the word indicating that the data set is likely to have a 
parameter of a value close to one being sought statistically. 

5 Swtch: "Switching parameter"; A two valued parameter whose purpose is to evaluate 
deRxInsl relative to some domain limits and shift between the two "domain" equations 
for deRxInsl(i). (Multiple Days' Version) 

Target: General: A number used in a feedback algorithm representing the desired result. 

10 

TargetAM: "Target Blood Glucose After Meals". Corrective insulin is calculated 
according to how high a patient's BG is above the TargetAM. See CorBol 

TargetTB: "Target Blood Glucose Before Meals". Corrective insulin is calculated 
15 according to how high a patient's BG is above the TargetTB. See CorBol 

Target Basal-to-Total Ratio: The desired value of daily Basal Insulin divided by Total 
Daily Insulin. Used in a feedback algorithm. 

20 TBCorIns(d-l,i+l): "Time-Boundary Corrective Insulin." The patient is usually 
instructed to test his or her BG at the time boundaries. This parameter is the previous 
day's total of Corrective Insulin taken on or near the time boundary, time(i+l). (Daily 
Update Version) 

25 TBCorIns(i+l): " Time-Boundary C orrective Insulin in the (i+l)th time interval." The 
patient is usually instructed to test his or her BG at the time boundaries. This parameter 
is the average over several days of the total of Corrective Insulin taken on or near a time 
boundary, time(i+l). (Multiple Past Days Version) 

30 TBCorTot(d): "Time-Boundary Corrective Insulin Total for the current day". (Daily 
Update Version) 

TBCorTot: "Time-Boundary Corrective Insulin Total": The sum of TBCorlns(i) over the 
whole day. (Multiple Past Days Version) 

35 

TDD(d-l): "Total Daily Dose" of insulin. The total amount of insulin a patient received 
during the previous day. (Daily Update Version) 

TDDavg: "Total Daily Dose of Insulin, average." The average of TDD for several days. 
40 (Multiple Past Days Version) 

Time(i): The "ith" time boundary 
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Time-Boundary Corrective Insulin: Corrective Insulin Boluses taken at or near a 
boundary of a time interval. Typically a BG test is taken at the time boundary (just 
before eating if the next interval is a meal interval) and the corrective bolus is calculated 
from this BG test. 

5 

TimeLabel(i): A short text phrase labeling each time boundary. PumpMaster offers 
several standard entries: Mid-Sleep, Pre-Breakfast, Pre-Lunch, Pre-Supper, Bedtime, 
Snack, and Basal Rate Change. 

10 Tmax: The last time boundary of a patient's day. Usually bedtime. 



4. GENERIC DERIVATIONS: 

15 4. 1 MAKING CHANGES IN THE CARBOHDRATE-TO-INSULIN RATIO, CIR: 

We can re-arrange the definition of CIR in two different ways, one for current conditions 
and one as a recommendation, assuming that the carbohydrates do not change: 

20 Carbs = CIR * Meal Insulin (21) 

CIRrec= Carbs / Recommended Meal Insulin (22) 

If Meal Insulin or Carbs are not known, then differential methods can be employed: 

CIRrec = CIR + dCIR (23) 
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30 
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We need the derivative of CIR with respect to Meal Insulin to use in the formulas below 
(the derivative is in parentheses): 



dCIR= (dCIR/dMeal Insulin) * dMeal Insulin (24) 
And 

35 dMeal Insulin = dCIR / (dCIR/dMeal Insulin ) (25) 

The following derivation leads to the formula for calculating dCIR, the change in CIR, 
that will produce a change d(Meal Insulin) in Meal Insulin. There are several ways to 
calculate or estimate the derivative. To obtain them, start with the definition of CIR: 



CIR = Carbs / (Meal Insulin) (26) 
Differentiate using calculus: 
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(dCIR/dMeal Insulin) = - Carbs / (Meal Insulin) 2 (27) 

Substitute equation (21) into equation (27) to obtain: 

(dCIR/clMeal Insulin) = - CIR / (Meal Insulin) (28) 

This is an estimate of the derivative that depends on CIR. It is possible to eliminate CIR 
from the equation by substituting a correlation from a statistical study. For example, the 
AIM correlation, (Ref 1) is used below. It says: 



CIR = Kcir * BodyWt / TDDavg (29) 
where Kcir is a statistically-determined constant. 
15 Substitute this into equation (28) for CIR to obtain the following generic equation: 

(dCIR/dMeal Insulin) = - ( Kcir * BodyWt / TDDavg ) / (Meal Insulin) (30) 
If Meal Insulin is unavailable, it can be calculated as: 

20 

Meal Insulin = TDDavg - Basal - Corrective Insulin (31) 
Leading to: 

25 (dCIR/dMeal Insulin) = - ( Kcir * BodyWt / TDDavg ) / (TDDavg - Basal - Corrective 
Insulin) (32) 

The AIM statistical studies show that Basal is slightly less than half of TDDavg. 
Corrective Insulin is usually small, so we can say that Meal Insulin is approximately half 
30 of TDDavg. This leads from equation (30) to the equation: 

(dCIR/dMeal Insulin) = - ( Kcir * BodyWt * 2 / TDDavg 2 ) (33) 

and finally: 

35 

dCIR = - ( Kcir * BodyWt * 2 / TDDavg 2 ) * dMeal Insulin (34) 
Another correlation (Ref 2) estimates CIR= Kcirw/TDD. By similar steps this leads to: 
40 dCIR =-( Kcirw * 2 / TDDavg 2 )* dMeal Insulin (35) 

The change in Meal Insulin, dMeallns, may be obtained as follows: 
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dMeallns = MeallnsNew - MeallnsOld (36) 

dMeallns = Carbs * (1/CIRnew - 1/CIRold ) (37) 

5 or it may be input by the Practitioner. 

4.2 SOME IMPORTANT EQUATIONS: 

4.2.1 CONVENTIONAL INSULIN NOMENCLATURE 
10 The earlier section showing the different types of insulin is re-written below using the 
variable names. 

Recall the definitions: Rxlnsl: Prescription Insulin, 

TBCorTot: Total Day's Time-Boundary Corrective Insulin, 
AMCorTot: Total Day's After-Meal Corrective Insulin 
15 CorTot: Totol Day's Corrective Insulin 

TDD = Rxlnsl + CorTot (38) 
Where Rxlnsl = Basal +Meallns (39) 

20 For small changes: 

dRxInsl = dBasal + dMeallns (40) 

Within a time interval, similar equations holds true: 

25 

Rxlnsl(i) = Basal(i) + Meallns(i) (41) 
dRxInsl(i) = dBasal(i) + dMeallns(i) (42) 

The present invention incorporates the following concept: The Total Day's Corrective 
30 Insulin, CorTot represents an "error" in the patient's Prescription Insulin; the goal is 
CorTot=0. The practitioner (or an automatic algorithm) decides how much of CorTot to 
eliminate. This amount is called dRxInsl mimicking the words "Change in Total 
Prescription Insulin". By making this change, the program can indirectly cause CorTot to 
be reduced in the days ahead. dRxInsl may be determined differently, depending on the 
35 Version: 

Daily Update Version: Krxlnsl is a pre-determined constant <=1 . It is used to set 
dRxInsl in the manner shown below: 

40 dMeallns + dBasal = dRxInsl = Krxlnsl * CorTot (43) 

Similarly, in the time interval: 

dMeallns(i) + dBasal(i) - dRxInsl(i) = Krxlnsl * CorIns(i+l) (44) 
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At the present time Krxlnsl = 0.16 for the Daily Update algorithms, but the value of 
Krxlnsl is subject to adjustment for optimum safety and performance. 



5 Multiple Days' Data Version: In the Manual Sub- Version, dRxInsl is input by the 
Practitioner for each Patient/Practitioner Interaction. Any value is allowed up to a 
maximum of KrxInslMax*CorTot, where KrxInslMax is a fractional constant: 



10 

Each time interval gets its share of the "fix" in proportion to its share of the total error : 



15 At the present time, KrxInslMax = 0.5 but is subject to adjustment for optimum safety 
and performance. 

In the automatic Sub-version, an automatic method determines dRxInsl (explained in the 
section entitled" Automatic Multiple Days' Data (Digital Advisor)". 

20 The left-hand-side of equation (44) or (46) is implemented for the interval in one of two 
ways. These are given the name "floats": 

The Meallns Float: dBasal(i) is estimated first; then the invention calculates 

dMeallns(i) = dRxInsl(i) - dBasal(i). Some of the ways for 
25 estimating dBasal(I) are: 



dMeallns + dBasal = dRxInsl <= KrxInslMax * CorTot 



(45) 



dMeallns(i) + dBasal(i) = dRxInsl(i) = dRxInsl * (CorIns(i+l) /CorTot) 



(46) 



Borrowing new Basal Rate from another interval; 
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Estimating Basal Rate from another algorithm (e.g. a Basal 
Float); 
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Estimating dBasal(i) from the carb-free latter part of a time 
interval using TBCorlns(i) as the error indicator. The float is 
done on the first part of the interval, using AMCorlns(i) as the 
error. (47) 



or 
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The Basal Float: 



dMeallns(i) is estimated first; then the invention calculates 
dBasal(i) = dRxInsl(i) - dMeallns(i) . Some of the ways for 
estimating dMeallns(i) are: 



Borrowing Meallns(i) or CIR(i) from another interval; 



Estimating dMeallns(i) as a share of day's total dMeallns in 
the same proportion as CarbSh(i)/CarbShTot; 

Estimating dMeallns(i) as a share of day's total dMeallns in 
5 the same proportion as MealIns(i)/MealInsTot (48) 

4.2.2 ENHANCED INSULIN NOMENCLATURE: 



The types of insulin-delivery systems that incorporate Insulin-on-board calculations make 
10 it safer for the patient to use After-Meal Corrective Insulin dosing, (variable name: 
AMCorlns) It i s c ailed " corrective" b ecause i t u ses t he correction f ormula, equation 
(14). However, since it is part of the insulin needed to cope with a specific meal of carbs, 
it is more convenient to lump it together with Meal Insulin by defining new variables: 



15 Enhanced Meal Insulin: eMeallns = Meallns + AMCorlns 

Enhanced Prescription Insulin: eRxInsl = Basal + eMeallns 

This is explained better as follows: 

TDD = eRxInsl + TBCorlns (49) 

20 

where eRxInsl = Basal + eMeallns (50) 

and eMeallns = Meallns + AMCorlns (5 1 ) 

25 The two quantities in eMeallns are added together at the start of these calculations and 
treated as a single variable. 

For small changes: 

30 deRxInsl = dBasal + deMeallns (52) 
Within a time interval, similar equations holds true: 
eRxInsl(i) = Basal(i) + eMeallns(i) (53) 

35 

deRxInsl(i) = dBasal(i) + deMeallns(i) (54) 



The present invention incorporates the following concept: The Total Time-Boundary 
Corrective Insulin, TBCorTot represents an "error" in the patient's Enhanced Prescription 
40 Insulin; the goal is TBCorTot=0. The practitioner (or an automatic algorithm) decides 
how much of TBCorTot to eliminate. This amount is called deRxInsl mimicking the 
words "Change in Total Enhanced Prescription Insulin". By making this change, the 
program can indirectly cause TBCorTot to be reduced in the days ahead. deRxInsl may 
be determined differently, depending on the Version: 
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Daily Update Version: Krxlnsl is a pre-determined constant <=1. It is used to set 
deRxInsl in the manner shown below: 

deMeallns + dBasal = deRxInsl = Krxlnsl * TBCorTot (55) 

Similarly, in the time interval: 

deMeallns(i) + dBasal(i) = deRxInsl(i) = Krxlnsl * TBCorIns(i+l) (56) 

At the present time Krxlnsl = 0.16 for the Daily Update algorithms, but the value of 
Krxlnsl is subject to adjustment for optimum safety and performance. 

Multiple Days' Data Version: In the Manual Sub- Version, deRxInsl is input by the 
15 Practitioner for each Patient/Practitioner Interaction. Any value is allowed up to a 
maximum of KrxInslMax*TBCorTot, where KrxInslMax is a fractional constant <=1. 

deMeallns + dBasal = deRxInsl <= KrxInslMax * TBCorTot (57) 

20 Each time interval gets its share of the "fix" in proportion to its error : 

deMeallns(i) + dBasal(i) = deRxInsl(i) = deRxInsl * (TBCorIns(i+l) /TBCorTot) (58) 

At the present time KrxInslMax = 0.5 but is subject to adjustment for optimum safety and 
25 performance. In the automatic Sub-version, an automatic method determines deRxInsl 
(explained "Automatic Digital Assistant). 

The left-hand-side of equation (56) or (58) is implemented for the interval in one of two 
ways. These are given the name "floats": 

30 

The Meallns Float: dBasal(i) is estimated first; then the invention calculates 

deMeallns(i) = deRxInsl(i) - dBasal(i). Some of the ways for 
estimating dBasal(i) are: 

35 Borrowing Basal(i) from another interval 

Estimating Basal(i) from another algorithm 
(e.g. a Basal Float) (59) 
or 

40 

The Basal Float: deMeallns is estimated first; then the invention calculates 

dBasal(i) = deRxInsl(i) - deMeallns(i) Some of the ways for 
estimating deMeallns are: 
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Borrowing eMeallns(i) or CIR(i) from another interval 

Estimating deMeallns(i) as a share of day's total deMeallns in 
the same proportion as CarbSh(i)/CarbShTot. (60) 

5 

Estimating deMeallns(i) as a share of day's total deMeallns in 
the same proportion as eMealIns(i)/eMealInsTot 

4.2.3 GENERAL 

10 

For a given Version, the different time intervals may have different types of float. 

The Basal float has the advantage that it can be used for all time intervals in the day, 
enabling easy changes to the patient's daily routine. 

15 

The Meal Insulin Float has the advantage that Basal Rate schedule can be kept simple, 
and changes in eating habits can be handled by CIR which addresses the changes only if 
the p erson consumes carbs and takes a meal bolus. However, since the Meal Insulin 
Float can only affect the t ime i ntervals w ith m eals in them, the Basal Float equations 
20 must still be used for the non-meal intervals. 

4.3 THE TWO-LEVEL BASAL SYSTEM, and the Kf CALCULATOR 

Some Practitioners prefer a "Two Level Basal" schedule, that assigns a single Basal Rate, 
25 BRf, (for BR fasting) around the clock except for the Late-Sleep time interval, just before 
waking. The Late-Sleep interval is different because of the "Pre-Dawn Phenomenon", 
which is the requirement for somewhat more insulin during the second half of the sleep 
period. This results in the two-levels of Basal Rate: BRf and BRlateSlp. Finding BRf 
becomes an important task. 

30 

A non-meal interval would be best for determining BRf. Typically, the non-meal 
intervals are the Early-Sleep and Late-Sleep intervals. H owever, neither of these is a 
perfect candidate: 

The Early-Sleep interval, if truly meal-free, would give the best value of BRf, but often 
35 patients "cheat" by taking a bedtime snack. Also, the patient must interrupt his or her 
sleep in order to take the necessary mid-sleep BG test, so the data is often not available. 
The Late-Sleep or pre-dawn interval is the most readily available and most often meal- 
free, but this interval is unsuitable because of the "dawn phenomenon". Neither of these 
is a perfect candidate for BRf. 

40 

The present invention incorporates two generalized methods that are designed to provide 
several ways to calculate BRf Some nomenclature should be explained: BRf is intended 
to be used in most intervals excluding a few which are called "excluded intervals". One 
interval or category of BR's with good statistical reliability may be nominated as a 
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"reliable" interval (with basal rate BRreliable). Another interval or data set may be 
nominated as a "similar" interval because its basal rate, BRsimilar, is similar to BRf. The 
present invention uses two generic methods to obtain BRf. 

5 4.3.1 Statistical Correlation: 

BRf may be obtained by a familiar least squares formula that statistically correlates 
paired values of "reliable" and "similar" data over the recent calendar period with a 
correlation constant, Kf, estimated to fit the formula BRsimilar = Kf * BRreliable. 

10 Kf = SUM[ ( BRsimilar ) * (BRreliable) ] / SUM[ BRreliable 2 ] (61) 

Once calculated, the correlation constant, Kf, is used to convert the values of the 
"reliable" basal rates into values of BRf until the next Practitioner interaction: 
BRf = BRreliable * Kf (62) 

15 

4.3.2 Averaging: 

BRf may be obtained by averaging certain sets of BRsimilar data from the same day or 
the previous day (in the case of the Daily Update Version), or from the recent calendar 
period (in the case of the Multiple Days' Data Version). The smoothing effect of the 
20 averaging process helps to produce good values. 

Some particular applications of these two methods are described below: 

The the Early-Sleep BR's of nights in which there was no bedtime snack are used as 
25 BRsimilar. Paired with these are values of Late Sleep BR immediately following, in 
the early morning of the next calendar day, which are used as BRreliable. The 
formula for subsequent use is: BRf = Kf * BRlateSlp (63) 

The average of BR's except BRlateSlp may be nominated as the similar data set, 
30 paired with BRlateSlp as the reliable data set. The resulting formula appears the same, 
as equation (63) though the Kf is different. 

The average of all BR's except BRlateSlp may be used directly as BRf. 

35 The average of all BR's i.e. (Basal/24) may be used directly as BRf. 

The average of all BR's i.e. (Basal/24) may be nominated as the similar data set, 
paired with BRlateSlp as the reliable data set. Once again the resulting formula is the 
same as equation (63). 

40 

The Early-Sleep BR's of nights in which there was no bedtime snack are used as 
BRsimilar. Paired with these are "reliable" values of Basal/24. The formula for 
subsequent use is: BRf - Kf * Basal/24 (64) 
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The average of BR's except BRlateSlp may be nominated as the similar data set, 
paired with the average of all BG's, i.e. Basal/24 as the reliable data set. The resulting 
formula is the same as (64) but Kf is different. 

5 The average of the Basal Float 1 results for all the other intervals (except Late-Sleep 
interval) may be nominated as the similar data set, paired with BRlateSlp as the 
reliable data set. Once again the resulting formula looks the same as equation (63) 

The average of the Basal Float 2 results for all the other intervals (except Late-Sleep 
10 interval) may be nominated as the similar data set, paired with BRlateSlp as the 
reliable data set. Once again the resulting formula looks the same as equation (63) 

The average of the Meal Insulin Float 2 results for all the other intervals (except Late- 
Sleep interval) may be nominated as the similar data set, paired with BRlateSlp as the 
15 reliable data set. Once again the resulting formula looks the same as equation (63) 

Additional applications of this type are mentioned throughout the text that follows. 

General discussion of Kf: 
20 The Kf Calculator is particularly well-suited for the Multiple Days' Data Versions so that 
it will be performed external to the insulin-delivery device. It is designed to be operated 
by the Practitioner at the time of a patient/practitioner interaction. This is because it 
depends upon digital memory of a sufficient number of nights to obtain enough data 
points to calculate an accurate value for Kf. 

25 

5. TYPES OF INSULIN DELIVERY SYSTEMS: 

The following is a list of insulin delivery systems showing the Versions of the present 
invention applicable to each: 

30 

Type A pumps: Have memory for TDD for several days. They have memory of the 
programmed Basal schedule, Basal(i). They have no memory for CF, CIR, CIR(i), 
BG(i), MealBol(t), Meallns(i), TBCorlns, AMCorlns, CorBol(t), or Corlns(i). Therefore, 
average combined boluses for the day must be obtained by subtracting Basal from 
35 TDDavg. Example: Medtronic MiniMed Paradigm 511. 

Version 6.2. 1 . 1 .4: which uses Multiple Days' Accumulated Data. Program 
is external to pump. 

Type B pumps: In addition to the Type A memories, these have BG(i), actual carbs for 
40 CarbSh(i); CIR(i), CF, and Combined Boluses, Boli(i). Examples: Medtronic MiniMed 
Paradigm 512 and 712. 

Version 6.2.1.1.3: which uses Multiple Days' Data. Program external to pump. 
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Type C pumps: In addition to other Type B memories, these have memories for 
Meallns(i) and Corlns(i). Example: Deltec Cozmo. 

Version 6.2. 1 . 1 .2: which uses Multiple Day's Data. Program external to pump. 

5 Type D pumps: In addition to Type C memories, these pumps have, as a minimum, 
memories for AMCorlns(i) and TBCorlns(i). 

Versions 6.2.1.1,1: which uses Multiple Days' Data. Program external to pump. 

Type E pumps: Have all the parameters of the Type D pumps as a minimum and have the 
10 program installed internally. 

Versions 6.1.2: which uses automatic Daily Update. Program internal to pump. 
Version 6.2.1.1.1: which uses Multiple Days' Accumulated Data. Program 
external to pump. 

15 Subcutaneous or Inhaled Insulin Delivery in which the insulin delivery device and the 
BG meter clip together into a kit or otherwise communicate with each other so that data is 
recorded digitally. 

Versions 6.1.3: which uses automatic Daily Update. Program internal to the kit. 
Version 6.2. 1.2: which uses Multiple Days' Accumulated Data. Program external 
20 to the kit. 

For the purpose of ease of presentation, these are re-organized in the Table of Contents. 
The different types of algorithm are shown as well: 

25 6. DESCRIPTION BY. VERSION. EMBODIMENT. ALGORITHM: 



6.1. VERSIONS USING DAILY UPDATE: 
30 6.1.1 PRELIMINARY DERIVATIONS: 

6.1.1.1 TIME INTERVALS AND CORRECTIVE INSULIN 
The index "d" denotes the present day. 

The index "i" denotes the "ith" time boundary or time interval (following the time 
35 boundary). 

The Time-Boundary Corrective Insulin, TBCorIns(d,i) is the sum of corrective boluses 
occurring at the beginning of the ith interval Each bolus needs to be identified as a 
Time-Boundary bolus, and it needs to be identified with the ith interval. This may be 
40 accomplished by a suitable combination of the following methods: 

identified with the ith interval if it falls between the midpoint of the previous 
interval and the midpoint of the ith interval. 
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identified as a Time-Boundary Corrective bolus if it comes before the Meal Bolus 
with the same index number. 

one or both of these identifications input by the patient at bolus-time, using the 
5 controls on the insulin-delivery device. 

The After-Meal Corrective insulin AMCorIns(d 5 i) is the sum of corrective boluses, 
identified by a suitable combination of the following methods: 

10 identified with the ith interval if it falls within the ith interval. 

identified as an After-Meal Corrective bolus by the fact of coming after the Meal 
Bolus for the meal that marks the start of the interval. 

15 one or both of these identifications input by the patient at bolus-time, using the 

controls on the insulin-delivery device. 

6.1.1.2 Krxlnsl: GOVERNS THE SIZE OF INSULIN CHANGES 

20 Suppose the operation of the invention is analyzed starting at an arbitrary starting time, 
when TBCorTot = TBCorTotStart. This represents an error to the invention. The 
purpose is to decrease each day's "error", TBCorTot, by a fractional amount 
KrxInsPTBCorTot, where KrxInsl<=L This means that, starting with a "starting error", 
TBCorTotStart, the error should decrease in a geometric sequence. After a number of 

25 days, Ncycle, 

The remaining error, TBCorTot(Ncycle) = TBCorTotStart* [ 1 -KrxInsl] Ncycle (65) 

This remaining error approaches near-zero as Ncycle increases. The phrase "near-zero" 
can be defined by setting a reasonably small number for a Percent Remaining Error. 
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Percent Remaining Error = [l-KrxInsl] Ncycle (66) 

For a suitable combination of Percent Remaining Error and Ncycle, the equation can be 
solved for Krxlnsl. This is done by the invention as follows: 

Krxlnsl = 1 - ( Percent Remaining Error f mc ^ (67) 



The Percent Remaining Error is set by the inventors but is subject to change. (A typical 
value is 10%). The choice of Ncycle is left to the practitioner. (A typical value is 14 
40 days). This example leads to an exemplary value of Krxlnsl=0.16 . In plain language, 
"A choice of 0.16 for Krxlnsl eliminates 90% of the error in two weeks." It is used in 
the Daily Update algorithms as follows: 

deRxInsl = Krxlnsl * TBCorTot (68) 
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or 

dRxInsl = Krxlnsl * TBCorTot (69) 
5 In the interval: 

deRxInsl(i) = Krxlnsl * TBCorIns(i+l) (70) 

or 

10 dRxInsl(i) = Krxlnsl * TBCorIns(i+l) (71) 
6, 1.2 FOR PUMPS (TYPE E) 
6.1.2.1 BASAL-TO-TOTAL RATIO: 

15 

The Basal(d)/TDD(d) ratio is called BoT(d) for "Basal over TDD". It is calculated every 
day. It has been determined by medical studies that certain ratios of Basal(d)/TDD(d) 
lead to better management of diabetes, a target for Basal(d)/TDD(d) ratio can be set into 
the p rogram b y t he p atient o r p ractitioner. It i s c ailed BoTTgt f or "Basal over TDD, 
20 Target". 

The invention contains a feedback factor to bring BoT(d) to BoTTgt. The feedback 
factor incorporates a constant Kfbk to regulate the speed of convergence: 

25 BoTFbk(d) = BoT(d) + Kfbk*( BoTTgt - BoT(d) )*sign(deRxInsl) (72) 

Where the resulting value is not allowed to be less than zero or greater than one. 

This factor, multiplied times the day's proposed total insulin change will yield the 
30 amount of change to Basal. The "sign" function ensures that the magnitude of the change 
is correct for the "direction" of the change. The constant Kfbk is adjusted to achieve the 
optimum speed of convergence to BoTTgt. The feedback factor is applied as follows: 
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dBasal - BoTFbk(d) * deRxInsl (73) 

deMeallns - (l-BoTFbk(d)) * deRxInsl (74) 

The invention provides a recommendation for BoTTgt, named BoTTgtRec, which is 
calculated as follows: 

BoTTgtRec = 1 - 4*CarbShTot*(Average Glycemic Index)/(a statistical correlation for 
caloric intake as a function of height, weight, and other easily known patient parameters) 
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This formula changes carbs to calories by use of the conversion factor 4, then multiplies 
by glycemic index to obtain the immediately-available calories from the carbs, then 
divides by the patient's caloric intake as estimated by body conformation. This gives the 
requirement for Meallns/TDD. Basal to Total ratio is one minus this quantity. 

5 

The present versions of the invention generally do not use this concept for both dBasal 
and deMeallns in the same time interval, but instead use one of them and determine the 
other parameter by subtraction from the total, deRxInsl in the manner described in the 
definition of "Float". 

10 

6.1.2.1 BASAL FLOAT L ( deMeallns is proportional to eMeallns or Carbs) 



This algorithm estimates Meal Insulin using the feedback factor, BoTFbk(d). The 
15 Enhanced insulin nomenclature system is used. The total change in Enhanced Meal 
Insulin is taken from equations (72) and (76): 

deMeallns(d) = Krxlnsl * (1 - BoTFbk(d-l)) * TBCorTot(d-l) (75) 

20 The portion of this quantity, which is assigned to the interval is proportional to the 

interval's share of Enhanced Meal Insulin from the day before. This concept is applied in 
the following equation: 

deMealIns(d,i) = Krxlnsl *(1- BoTFbk(d-l)) * TBCorTot(d-l)* eMealIns(d- 
25 1 ,i)/eMealInsTot(d- 1 ) (76) 

eMealIns(d,i) = eMealIns(d-l,i) + Krxlnsl *(1- BoTFbk(d-l)) * TBCorTot(d- 

1 )*eMealIns(d- 1 ,i)/eMealInsTot(d- 1 ) (77) 

30 Using equation (22) gives: 

CIR(d,i) = CarbSh(d-l,i) / [eMeallns(d-U) + Krxlnsl *(1- BoTFbk(d-l)) * TBCorTot(d- 
l)*eMealIns(d-l,i) / eMeallnsTot(d-l)] (78) 

35 As in the definition of a Basal Float, statement (60): 

BR(d,i) = BR(d-l,i) + (deRxInsl(d,i) - deMealIns(d,i))/dt(i) (79) 

Using equation (70) gives: 

40 

BR(d,i) = BR(d-l,i) + Krxlnsl * [ TBCorIns(d-l,i+l) - (l-BoTFbk(d-l))*TBCorTot(d- 
l)*eMealIns(d-l,i) / eMeallnsTot(d-l) ] / dt(i) (80) 
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This algorithm can be used in all time intervals "across the board". If used in all 
intervals, this algorithm has the effect of maintaining the original prescribed "shape" of 
the CIR or Meallns schedule; that is, each eMealIns(d,i) is multiplied by the same factor, 
so that they rise or fall in unison. 

5 

6.1.2.2 BASAL FLOAT 2, ( deMeallns value from outside the interval ) 

This algorithm is used i n i ntervals w here t he M eallns o r C IR v alue i s b orrowed from 
outside the interval (pegged to another interval.) This method is useful in non-meal time 
intervals, because meal insulin data in these intervals is not dependable enough for 
10 calculations. All the CIR's go up or down in proportion to the "Key CIR". The CIR 
schedule maintains the same "shape" and each CIR preserves the same ratio to the key 
CIR that it had on the first day, d=l. The key CIR may be chosen from several sources, 
usually from a time interval that has meals and uses a Meal Float algorithm. The number 
of the key CIR interval is keyc. The Enhanced insulin nomenclature system is used. 

15 

CIR(d,i) = [CIR(d-l, keyc) / CIR(1, keyc) ]* CIR(1, i) (81) 

Equation (37) is adapted as shown below. The Enhanced nomenclature system is used, 
but note that the desired change is the same in both systems: 

20 

deMealIns(d,i)= dMealIns(d,i) = CarbSh(d-l,i)*[l/CIR(d,i)- l/CIR(d-l,i)] (82) 
BR(d,i) = BR(d-l,i) + dBR(d,i) (83) 
25 The "Float" equation can be obtained from statement (60) as follows: 

BR(d,i) = BR(d-l,i) + (deRxInsl(d,i) - deMealIns(d,i))/dt(i) (84) 
Then equation (56) is applied as follows: 

30 

BR(d,i) = BR(d-l,i) + [ KrxInsl*TBCorIns(d-l,i+l) - deMealIns(d,i) ] / dt(i) (85) 

BR(d,i) = BR(d-l,i) +[ Krxlnsl * TBCorIns(d-l,i+l) - CarbSh(d-l,i)*[l/CIR(d,i)- 
l/CIR(d-l,i)]]/dt(i) (86) 

35 

6.1.2.3 MEAL INSULIN FLOAT 1 

This algorithm pegs BR to a "Key" basal rate, i.e. all the basal rates go up or down in 
proportion to the key basal rate. The B asal schedule maintains the same "shape" and 
40 each Basal Rate preserves the same ratio to the key basal rate that it had on the first day, 
d=l . The key basal rate may be chosen from several sources, usually from a time interval 
that has no meals and uses a Basal Float algorithm. The key interval number is keyb. The 
Enhanced insulin nomenclature system is used. 
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BR(d,i) = [BR(d-l, keyb) / BR(1, keyb) ]* BR(1, i) 



The key basal rate can also be obtained from a calculation that uses the overall result of a 
Basal Float 1 algorithm, i.e. Basal(d-l) plus the correction. The formula uses d=l like 
the one above. 

BR(d,i) = [Basal(d-l) + BoTFbk(d-l)*KrxInsl*TBCorTot(d-l)] / Basal(l) * 

BR(U) (88) 



The key basal rate may also be obtained from: 

10 

From a sum of basal rates over all the intervals. 

Basal(d) = SUM( BR(d,i)*dt(i) ) (89) 
BR(d,i) = [(Basal(d)/Basal(l)] * BR(l,i) (90) 

15 Whatever the source of BR(d,i), the change in Basal Rate is: 

dBR(d,i)= BR(d,i)-BR(d-l,i) (91) 

As in the definition of a Meal Insulin Float, statement (59): 

deMealIns(d,i) = Krxlnsl * TBCorIns(d-l,i+l) - dBR(d.i) *dt(i) (92) 

20 

deMealIns(d,i)=KrxInsl*TBCorIns(d- 1 ,i+l )-( BR(d,i) - BR(d- 1 ,i) )*dt(i) (93) 

eMealIns(d,i) = eMealIns(d-l,i) + Krxlnsl * TBCorIns(d-l,i+l) - (BR(d,i) - BR(d-l,i) ) * 
dt(i) (94) 

25 

Paraphrasing equation (22) gives: 

CIR(d,i) = CarbSh(d-l, i) / eMealIns(d, i) (95) 

CIR(d,i) = CarbSh(d-l,i) / [ eMealIns(d-l,i) + Krxlnsl * TBCorIns(d-l,i+l) - (BR(d,i) - 
30 BR(d-l,i))*dt(i)] (96) 



6.1.2.4 OVERVIEW OF 6. 1 . 1 .2 and 6. 1 . 1 .3 



35 The Basal Float 2 algorithm is used in conjunction with the Meal Insulin 1 float as 
follows: 

The practitioner's latest prescription is input for the first day, d=l. This data is the 
starting point for the automatic daily update sequence. The non-meal intervals are 
provided with Basal Float 2 algorithms, pegging their CIR's to a "Keyc" interval. One of 
40 the non-meal-containing intervals is nominated as the "Keyb" Basal Rate (usually the 
Late-Sleep interval or the average Basal/24). The Basal rates in the meal-containing 
intervals are pegged to the Key B asal Rate to maintain their original ratios to the key 
basal rate. The meal-containing intervals are provided with Meal Insulin Float 
algorithms. 



36 



6.1.2.5 MEAL FLOAT 2 (uses AMCorlns(i) as error ) 

This algorithm makes use of after-meal BG testing. These tests and a ssociated After- 
Meal Corrective Boluses divide the interval into two parts. The algorithm uses 
5 AMCorIns(d-l,i) as an "error" indicator for the first part of the interval, where the float is 
calculated. The value of dBasal(d-l,i) is obtained from the carb-free second part of the 
interval, which uses TBCorIns(d-l,i+l) as an error term for dBasal(d-l,i) The 
"Enhanced" insulin terminology is not used. 

10 Define: Timeb(d-l,i)= The time of the after-meal bolus after the ith time boundary. 

dtb(d-l,i) = Time(i+1) - Timeb(d-l,i) (97) 
dta(d-l,i) = Timeb(d-l,i) - Time(i) (98) 

15 In the second part of the interval: 

BR(d,i) = BR(d-l,i) + Krxlnsl * TBCorIns(d-l,i+l) /dtb(d-l,i) (99) 

Having obtained Basal Rate from the second part of the interval (or another source), the 
following statement can be said of the first part of the interval: 

20 

dBR(d,i) = BR(d,i) - BR(d-l,i) (100) 
Meallns(d,i) = Meallns(d-l,i) + Krxlnsl *AMCorIns(d-l,i) - dBR(d,i) * dta(d-l,i) (101) 
25 CIR(d,i) = CarbSh(d-l,i) / Meallns(d,i) (102) 



6.1.3 FOR MULTIPLE DOSE INJECTION (MDI) AND INHALED INSULIN: 
The Basal insulin is administered in the form of long-acting insulin as infrequently as 
30 once per day. Corrective Insulin and Meal Insulin are administered as needed in the form 
of injected or inhaled short-acting insulin. This algorithm is suitable for "kits" in which 
the BG Meter and the insulin delivery device clip together or are otherwise linked, in 
order to maintain a combined digital history and a place for the program to reside. 

35 6.1.3.1 BASAL FLOAT 1 . (not used with MDI or Inhaled Insulin) 

6.1.3.2 BASAL FLOAT 2, ( The value of dMeallns is from outside the interval ) 

This algorithm is used in intervals where the Meallns or CIR value is borrowed from 
40 outside the interval (pegged to another interval.) This method is useful in non-meal time 
intervals, because meal insulin data in these intervals is not dependable enough for 
calculations. All the CIR's rates go up or down in proportion to the "Key CIR". The 
CIR schedule maintains the same "shape" and each CIR preserves the same ratio to the 
key CIR that it had on the day of the most recent interaction with the Practitioner, d=l . 
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The key CIR may be chosen from several sources, usually from a time interval that has 
meals and uses a Meal Float algorithm. The number of the key CER interval is keyc. The 
Enhanced insulin nomenclature system is used. 



5 CIR(d,i) = [CIR(d-l, keyc) / CIR(1, keyc) ]* CER(1, i) (103) 

Equation (37) is adapted as shown below. The Enhanced nomenclature system is used, 
but note that the desired change is the same in both systems: 



10 deMealIns(d,i)=dMealIns(d,i) = CarbSh(d-l,i)*[l/CIR(d,i)- l/CIR(d-l,i)] (104) 

BR(d,i) = Basal(d-l)/24 + dBR(i) (105) 

from equation (54) it can be seen that: 
15 BR(d,i) = Basal(d-l)/24 + (deRxInsl(d,i) - deMealIns(d,i))/dt(i) (106) 

BR(d,i) = Basal(d-l)/24 + [ KrxInsl*TBCorIns(d-l,i+l) - deMealIns(d,i) ] / dt(i) (107) 

A non-meal interval using this algorithm (or some other suitable source) is used as the 
20 "reliable" data set. The resulting basal rate by equation (117) is used to determine the 
whole day's Basal. 

Basal(d) = 24 * (some conversion factor) * BR(d-l, reliable) (108) 

25 A recommended practice is the use of the Late-Sleep interval as the "reliable" interval. 
The f ollowing equation i s g iven w ithout e xplanation, a s i t w ill b e c overed b y another 
section. 



Basal(d) - 24 * Kf * BR(d-l,LateSleep) (109) 

30 

6.1.3.3 MEAL INSULIN FLOAT 1 

The Enhanced insulin nomenclature system is used. This algorithm is used primarily for 
meal-containing intervals. The whole day's Basal(d) can be from a key interval's Basal 
35 Float calculation as in equation (119). 



dBasal(d,i) = (Basal(d) - Basal(d-l)) * dt(i) / 24 (110) 
Whatever the source of Basal data, the Meal Insulin Float continues as follows: 

40 

The present invention calculates a schedule of new CIR(d,i) values, using the equation 
below, which is similar to equation (103): 
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CIR(d,i) = CarbSh(d-l,i) / [ eMealIns(d-l,i) + Krxlnsl * TBCorIns(d-l,i+l) - (Basal(d) - 
Basal(d-l))*dt(i)/24 ] (111) 



5 6.1.3.4 OVERVIEW OF 6.1.3 

The day's schedule is programmed as follows: An interval ( usually a non-meal interval) 
is nominated as the "reliable" interval and is provided with a Basal Float 2 algorithm, 
with its CIR pegged to a meal interval. The BR from this "reliable" interval is used to 
determine the round-the-clock BR and the total daily Basal dose, which equals a constant 

10 (e.g. Kf or a constant of that type) times BRreliable times 24. The meal-containing 
intervals are provided with Meal Insulin Float algorithms and use the same single-valued 
round-the-clock BR. The non-meal intervals other than the "reliable" interval contain no 
calculations; their CIR's are pegged to a meal interval and their BR's are the same as all 
the others. 



15 



25 



6.1.3.5 MEAL INSULIN FLOAT 2, (uses AMCorlns(i) as an error term) 



This algorithm makes use of after-meal BG testing. These tests and associated bolus 
divide the interval into two parts. The algorithm uses AMCorIns(d-l,i) as an "error" 
20 indicator for the first part of the interval, where the float is done. 

Define: Timeb(d-l,i)= The time of the after-meal bolus after the ith time boundary. 

dtb(d-l,i) = Time(i+1) - Timeb(d-l,i) (112) 
dta(d-U) = Timeb(d-U) - Time(i) (113) 



The Basal Rates are the same in all intervals; only Basal(d) is needed. It may be obtained 
from one of several sources. 



The carb-free second part of the interval uses TBCorIns(d-l,i+l) as an error term 
30 for dBasal(d,i) The "Enhanced" insulin terminology is not used. In the second 

part of the interval: 

BR(d,i) = BR(d-l,i) + Krxlnsl * TBCorIns(d-l,i+l) /dtb(d-l,i) (114) 

35 This result, obtained from a single "key" interval yields Basal(d): 

Basal(d) = [BR(d, key) / BR(1, key) ]* Basal(l) (115) 

The full schedule of results from the second part of the interval, equation (116) 
40 may be converted to a full day's basal as follows: 

Basal(d) = SUM over i [ BR(d,i) * dt(i) ] (117) 
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A Basal Float 1 algorithm like version 6.1.2.1 in a "reliable" interval may be 
used to provide a round-the-clock basal rate. 

A Basal Float 1 algorithm used in an overall manner, to obtain the day's total 
5 basal: 

Basal(d) = [Basal(d-l) + BoTFbk(d-l)*KrxInsl*TBCorTot(d-l)] (118) 

Whatever the source for Basal(d), the following statement can be said of the first part of 
10 the interval: 

Meallns(d,i) = Meallns(d-l,i) + Kxxlnsl *AMCorIns(d-l,i) - (Basal(d) - Basal(d-l) * 
dta(i)/24 (119) 

15 If the source is the second part of the interval, then this becomes: 

Meallns(d,i) = Meallns(d-l,i) + Krxlnsl *AMCorIns(d-l,i) - BoTFbk(d-l) 

* KrxInsl*TBCorTot(d-l)* dta(i)/24 (120) 

20 CIR(d,i) = CarbSh(d-l,i) / Meallns(d,i) (121) 

6.1.4 SKIPPED BG'S 

Versions 1 and all sub-versions keep track of a special "virtual interval" consisting of the 
25 combined Early-Sleep and Late-Sleep intervals. This is for use in the event that the 
patient skips the Mid-Sleep BG test. A Basal Float calculation is kept current over this 
combined interval and is used to provide the basal rate for both component invervals if 
the BG test between them is skipped. This technique is used for often-skipped BG's like 
the Mid-Sleep BG just described. This technique can be used for any number of adjacent 
30 time intervals with skipped BG's on the boundaries. However, for seldom-skipped BG's 
it is more expedient to merely substitute BG=Target so that a zero value is calculated for 
Corlns(i). This causes no change and allows the previous value to remain. F or this 
reason, the BG(d,i) parameters have default values of TargetAM or TargetTB until input 
is made. 

35 

6.1.5 A MODIFIED BOLUS CALCULATOR 

Many pump models provide a bolus calculator programmed inside. It performs the 
calculations in equations (20) and (14). The present invention treats exercise as a special 
40 adaptation of the Bolus calculator. There are input boxes for exercise in units of carbs, 
defining a new variable, ExerCarbs. The arithmetic difference between them is used as a 
variable Carbm shown below: 

Carbm = Carbs - ExerCarbs (122) 
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This may also be used for exercise alone without carbs. The modified Bolus Calculator 
also c an c alculate C orrection Boluses. They must be designated "Time-Boundary" or 
"After-Meal" by the patient to flag the memory record and to select the correct one of the 
5 two targets: Target AM and TargetTB. The modified boluses (meal and correction) are 
summed. If the result is positive, the pump infuses the calculated insulin amount as an 
ordinary bolus. If the result is negative, the pump suspends the Basal pumping for an 
amount of time calculated as follows: 

10 TimeOut= - NegativeBolus(t) / BR(d,i) (133) 

Or it may reduce the basal for a time: 

TimeReduced = - NegativeBolus(t) / (BR(d,i) - BRreduced ) (134) 

15 

6.1.6 CHANGING THE PATIENT'S SCHEDULE 

Section 1.1 above has described a normal 9 to 5 day. However, patients have many 
different schedules. To allow for this, the interval type is identified for each time interval 
20 by a parameter IntrvlType, which can have the following values: 

Mfor "Meal Interval" 

Kb for "keyb interval", typically the Late-Sleep interval, or other Non-Meal 
Interval in which food is never consumed. 
25 Kc for "keyc interval" the source of pegged CIR's 

S for "snack interval" in which food is occasionally consumed. 
R for "reliable" 

The adjustments to this feature are made by the Practitioner (not the patient unless it is 
with Practitioner guidance). The practitioner makes sure there is always a value of 
30 IntrvlType for each time interval. He or she can change this if desired. There is no 
restriction on the others. The present invention processes this data as follows: 

If IntrvlType=Kb, go to Algorithm 6.1.2.3, Basal Float 2 
This results in a calculation of BR(d,keyb). 

35 

If IntrvlType=Kc, go to Algorithm 6.1.2.4, Meal Insulin Float 4 
This results in a calculation of BR(d,keyb). 

If IntrvlType=R, go to Algorithm 6.1.2.3, Basal Float 2 
40 This results in a calculation of BRsimilar. 

If IntrvlType = S go to Algorithm 6.1.2.3, Basal Float 2 

IfIntrvlType = MTHEN: 
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IF a Time-Boundary BG is missing, go to Algorithm 6.1,4 

IF an After-Meal Corrective Dose is present, go to Algorithm 6.1.2.6 Meal, Insulin 
Float 2 

ELSE go to Algorithm 6.1.2.4, Meal Insulin Float 1 

5 

6.2 VERSIONS USING MULTIPLE-DAYS' DATA: 

Instead of re-adjusting the patient's parameters with the previous day's data, these 
Versions use the average accumulated data over the calendar period prior to the 

10 patient/practitioner interaction. The data for each time interval during the day is averaged 
separately over all the days. The day index, "d", is dropped and instead, the new 
parameters are distinguished from the current parameters as follows: Current parameters 
have no suffix. Calculated and recommended parameters have a suffix "rec". Prescribed 
parameters have suffix "rx". The "rx" parameters are input by the practitioner after 

15 considering the recommendations of the "rec" parameters. Parameters marked 
"Parameter(i)" are for the "ith" interval. Parameters with nothing in parentheses are non- 
scheduled parameters for which there is only one value for each patient/practitioner 
interaction. Thus "BRrec(i)" refers to the Recommended Basal Rate for the "ith" 
interval e.g. BRrecl, BRrec2 . . etc, and BRrx(i) refers to the Prescribed Basal Rates. 

20 

6.2.1 MULTIPLE DAYS' DATA (DIGITAL ADVISOR) FOR PRACTITIONERS 

Krxlnsl is replaced by deRxInsl/TBCorTot so (135) 
deRxInsl(i) - deRxInsl/TBCorTot * TBCorIns(i+l) (136) 

25 

Interactive Input Forms: 

An exemplary form of the present invention is a two-table Access Database. There is a 
"one-to-many" relationship between the "patients' table", Tp and the "interactions table", 
Ti. There is a digital interactive Input Form (described in more detail in section 6.2.1.1.1. 

30 

The practitioner obtains some of the patients data by the exam or interview process. 
Other data is obtained by downloading the data from the pump andBG meter, either 
directly into the invention or by using the manufacturer's software to make printed 
copies, which are then transcribed into the invention's interactive digital input form 
35 manually. The manufacturer's downloading software normally p rovides a schedule of 
Bins (time periods enveloping the primary Time Boundaries) and calculates averages 
pertaining to each time boundary for several parameters including BG(i), AMCorlns(i), 
TBCorlns(i) and Meallns(i). (Recall as an example, that Meallns(i) represents the 
average of total Meal Insulin boluses within the ith Bin.) 

40 

6.2.1.1 FOR PUMPS 

6.2.1.1.1 PUMP TYPE D and E 
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The formulas for pumps Types D and E are the simplest, so they will be discussed first. 
Interactive Input Forms: 

The main Input Form is the outer panel in Figure 1. It collects the patient's permanent 
5 demographic data for table Tp. A SubForm (inner panel) collects the data from the 
patient/practitioner i nteraction f or t able Ti. The Sub Form has two pages that can be 
reached by the scroll bar. The first page (Figure 1) addresses un-indexed parameters 
which have a single value for each patient/practitioner interaction. The second page 
(Figure 2) addresses the standard modal day's schedule, containing the parameters with 
10 time-interval indexes like those referred-to herein in the manner of "parameter(i)" 



6.2. 1 . 1 . 1 . 1 Basal Float 1 , (deMeallns is proportional to eMeallns or Carbs) 
This algorithm is similar to Basal Float algorithm 6.1 .2.2 for Daily Update Version. The 
15 Enhanced insulin nomenclature system is used. 

Transpose equation (52) to show that dBasal is determined by the other two parameters: 
dBasal = deRxInsl - deMeallns (137) 
dBasal(i) ^ deRxInsl(i) - deMeallns(i) (138) 

20 

deMeallns(i) is estimated by saying that the distribution of deMeallns among the 
intervals is proportional to the distribution of eMeallns, i.e.: 

deMeallns(i) = deMeallns * eMealIns(i)/eMealInsTot (139) 

25 

eMeallnsRec(i) = eMeallns(i) + deMeallns *eMealIns(i)/eMealInsTot (140) 
By equation (22): 

30 CIRrec(i)= CarbSh(i)/eMealInsRec(i) (141) 

ClRrec(i) = CarbSh(i)/[ eMeallns(i) + deMeallns *eMealIns(i)/eMealInsTot] (142) 

Note that both eMeallnsRec(i) and CIRrec(i) are merely the original value times an 
35 across-the-board factor. In the manner of statement (60), Basal is "Floated" by 
subtracting deMeallns(i) from deRxInsl(i): 

BRrec(i) = BR(i) + [ deRxInsl(i) - deMeallns(i)] / dt(i) (143) 

40 Substituting equations (133) and (136) into this gives: 

BRrec(i) = BR(i) + [ deRxInsl/TBCorTot * TBCorIns(i+l) 

- deMeallns* eMealIns(i)/eMealInsTot] / dt(i) (144) 
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Also calculated is the value of: BRaveRec = (Basal + deRxInsl - deMeallns ) / 24 

(145) 

The practitioner uses this version of the invention as follows: 
5 When all the patient's data has been input from downloads or otherwise, The Basal Float 
is almost ready to calculate the main goals (a schedule of recommended basal rates, 
BRrec(i) and a schedule of recommended CIRrec(i)). The practitioner needs to input at 
least two of the three quantities in equation (134). The present invention will do the rest. 
So, the practitioner looks over the data obtained so far, particularly TBCorTot. Then he 

10 or she makes a judgment as to "How much of TBCorTot do I want to add to Prescription 
Insulin as a change?" Then he or she inputs deRxInsl, which must be within the built-in 
input limits (see section: Limited Domain). Basal(i) will be "floated", so the other 
quantity needed from the practitioner is deMeallns. Before inputting he must ask 
himself, "When I prescribe this change, deRxInsl, how much of it do I want to assign to 

15 Enhanced Meal Insulin?" He can make this judgement by comparing the ratio of Basal / 
TDDavg (known as BoT in the present embodiment) to the optimum value from the AIM 
statistical studies (48% from latest publication). For instance, if Basal is too high he can 
use greater than half of deRxInsl as deMeallns. This will raise deMeallns relative to 
deRxInsl, thus lowering Basal. There may be other considerations of a medical nature 

20 that may influence the practitioner's decision. The choice of deMeallns may have input 
limits. The outputs are BRrec(i), BRavgRec, and CIRrec(i) as calculated by the 
equations above. The practitioner considers these recommended values and inputs the 
"rx" values based upon his judgment. 

25 

6.2.1.1.1.2 Basal Float 2, (not used with Multiple Days' Data) 

6.2.1.1.1.3 Meal Insulin Float 1 

30 This Version is similar to the Version 6.1.2.4, Meal Insulin Float. For input, it requires 
the Prescribed Basal Rates, BRrx(i), which are input by the Practitioner. The Practitioner 
may desire "advice" before inputting BRrx(i). The values of BRrec(i) and BRaveRec 
from the result of Version 6.2.1.1.1.1, Basal Float 1, above, are good advice, so they are 
provided on the same computer screen Input Form. The value of BRf from the Kf 

35 calculator is also good advice. The Enhanced insulin nomenclature system is used. 

Transpose equations (52) and (54) to show that deMeallns is determined by the other two 
parameters: 

deMeallns = deRxInsl - dBasal (146) 
40 deMeallns(i) = deRxInsl(i) - dBasal® (147) 

BRrx(i) = input by practitioner (148) 

The "Float" is very similar to equation (100); it is set up as follows: 



44 



Equation (132) is used to modify equation (100) with the result below: 
deMeallns(i) = (deRxInsl/TBCorTot)* TBCorIns(i+l) - (BRrx(i) - BR(i) )*dt(i) 

(149) 

5 

eMeallnsRec(i) = eMeallns(i) + deRxInsl/TBCorTot * TBCorIns(i+l) - (BRrx(i)-BR(i) ) 
*dt(i) (150) 
Using equation (22) gives: 

10 CIRrec(i) = CarbSh(i) / eMeallnsRec(i) 

(151) 

CIRrec(i) = CarbSh(i) / [eMeallns(i) + deRxInsl/TBCorTot * TBCorIns(i+l) - (BRrx(i) 
-BR(i))*dt(i)] (152) 

15 

The Practitioner considers this schedule of CIRrec(i) and then inputs: 

CIRrx(i) = input by the Practitioner (153) 

20 6.2. 1 . 1 . 1 .4 Overview of 6.2. 1 . 1 . 1 . 1 and 6.2. 1 . 1 . 1 .3 Multiple Days' Data, Pump Type D 

Recommendations are provided for Basal Float and Meal Insulin Float on the same Input 
Form (see Figure 2), so that the practitioner can consider Basal Float recommendations, 
BRrec(i), BRaveRec, and BRf when filling in the prescribed BRrx(i) schedule. The Meal 
25 Insulin Float uses these "BRrx(i)" values as input. Then it calculates CIRrec(i) as output 
for the meal-containing intervals. Then the Practitioner fills in the blanks for CIRrx(i). 
Backtracking the discussion a little: The practitioner's choice of BRrx(i) is very 
influential. Discussion: 



30 The practitioner can use all of the BRrec(i) values as BRrx(i) if desired. This has 

the effect of changing the entire meal insulin (or CIR) schedule by the same factor 
across the board. 

or 



35 A simple basal schedule can be used to give the patient more meal-skipping 

flexibility, as mentioned earlier. For instance: 



The practitioner can use the BR's for non-meal intervals "as is". This 
probably will include the "Early-Sleep" and "Late-Sleep" intervals. 

The Practitioner can use a single carefully-judged Fasting Basal Rate, BRf, to 
underlie the meal intervals. Determining this basal rate is one of the 
Practitioner's major tasks. The BRf Calculator may be used at the 
Practitioner's discretion. 
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10 



15 



6.2.1.1.1.5 Meal Insulin Float 2 (uses AMCorlns as the error ) 

This algorithm makes use of after-meal BG testing. These tests divide the interval into 
two parts. The algorithm uses AMCorlns(i) as an "error" indicator for the first part of 
the interval, where the float is done. The value of dBasal(i) is obtained from the carb-free 
second part of the interval, which uses TBCorIns(i+l) as an error term for dBasal(i). 
The "Enhanced" insulin terminology is NOT used. 

Define: Timeb(i)= The time of the after-meal bolus after the ith time boundary. 

dta(i) = Timeb(i) - Time(i) (154) 
dtb(i) = Time(i+1) - Timeb(i) (155) 

In the second part of the interval: 

BRrec(i) = BR(i) + Krxlnsl * TBCorIns(i+l) /dtb(i) (156) 

dBR(I) = BRrec(i) - BR(i) (157) 

20 Having obtained Basal Rate from the second part of the interval (or other source), the 
following statement can be said of the first part of the interval: 

MeallnsRec(i) = Meallns(i) + Krxlnsl *AMCorIns(i) - dBR(i) * dta(i) (158) 

25 If the second part of the interval was used, this becomes: 

MeallnsRec(i) = Meallns(i) + Krxlnsl *AMCorIns(i) 

- dta(i) * Krxlnsl * TBCorIns(i+l)/ dtb(i) (159) 

30 CIRrec(i) = CarbSh(i) / MealInsRec(d,i) ( 1 60) 

If desired, the factor Krxlnsl may be obtained from dRxInsl input by the practitioner as 
described earlier: 

35 Krxlnsl - dRxInsl / (AMCorTot+TBCorTot) (161) 

Alternatively, Krxlnsl may be replaced by two "K- factors", one for each part of the time 
interval: KrxInsAM and KrxInsTB: 

First, Basal Rate is determined from the second part of the interval: 
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BRrec(i) = BR(i) + KrxInsTB * TBCorIns(i+l) /dtb(i) (162) 
Then Meal Insulin is determined from the first part of the interval: 
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MeallnsRec(i) = Meallns(i) + KrxInsAM * AMCorlns(i) - dta(i) * KrxInsTB* 
TBCorIns(i+l)/ dtb(I) (163) 



CIRrec(i) = CarbSh(i) / MealInsRec(d,i) ( 1 64) 

5 

This alternative method may shift the Basal/TDD ratio and may be useful if such a result 
is intended. 

6.2.1.1.2 PUMP TYPEC 

10 Pumps of these types are intermediate between Types A and D. The algorithms of 
6.2.1.1.1.1 and 6.2.1.1.1.3 and their overview are applicable to TypeB pumps, except 
that: 

TBCorlns(i) is not available, so it must be calculated by the formula: 

15 

TBCorlns(i) = AVG over calendar period of ( (BG(t) - TargetTB)/CF) (165) 
AMCorlns(i) is not available so Meal Float 2 cannot be used. 

6.2.1.1.3 PUMP TYPEB 

20 Pumps of these types are intermediate between Types A and C. The algorithms of 
6.2.1.1.1.1, 6.2.1.1.1.3, and their Overview are applicable to Type B pumps, except that 
in addition to the limitations of Pump Type C: 

Meallns(i) is not available, so dMeallns(i) is estimated by: 

25 

dMeallns(i) = dMeallns * CarbSh(i)/CarbShTot (1 66) 

6.2.1.1.4 PUMP TYPE A 

Type A pumps are the simplest, but have complicated formulas for the reason that the 
30 values for the absent parameters must be calculated using estimation formulas. These 
estimation formulas add complexity. In addition to the limitations of Type B pumps, the 
Type A pumps are limited as follows: 

Because of the lack of data, the only appropriate algorithm for Type A pumps is 
35 "Basal Float 1 , similar to Version 6.2.1.1.1.1. 

There are no memories for BG(i), so values for TBCorlns(i) and AMCorlns(i) must 
be calculated from BG data downloaded from the BG meter. This data may contain 
anomalous or unused BG test results and is not as dependable as the BG data from the 
40 pump's memory. 

There are no memories for CarbSh(i), so the practitioner must estimate. The only use 
of the data is in the ratio CarbSh(i)/CarbShTot, so the units do not matter; the 
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practitioner may use grams, exchanges, percent of total, units of insulin, or any other 
units proportional to Meal Insulin. 

The lack of good data makes CIR hard to calculate. A single CIR is calculated for 
5 use round-the-clock. The differential formula (34) is used: 

dCIRcalcA = [ - ( Kcir * Wt *2 / TDDavg 2 ] * deMeallns (167) 

A slight complication is that Type A pumps have no memory for CIR, so the patient 
10 must keep track of it. To make this easier for the patient, only integer (or half integer 
values for CIR < 8) are prescribed. Therefore, a two-step manual input is used: 

dCIRcalcA is calculated by the invention and appears on the screen. The practitioner 
rounds it to the preferred rounded quantity and puts the value back into the reverse 
15 equation as dCIR. Then the invention calculates: 

deMeallnsCalcA = - dCIR / ( Kcir * Wt * 2 / TDDavg 2 ) (168) 

Then the practitioner inputs this value into the box for deMeallns. Not surprisingly, 
20 the resulting value of dCIRcalcA is integer-valued. The recommended value, 
CIRrecA, is automatically calculated as follows: 

CIRrecA = dCIRcalcA + CIR (1 69) 

25 

6.2.1.2 SUBCUTANEOUS MULTIPLE DOSE and INHALED INSULIN 



Manual injection patients usually use two types of insulin: 
30 Long-acting insulin for "basal" injections as infrequently as once per day 

Short-acting insulin for meal and correction boluses. 

The insulin delivery devices include "pens" and inhalers. 

35 There are "kits' currently being developed that consist of a BG meter that clips or links to 
an insulin-delivery device (an insulin injection "pen" or an insulin inhaler) in such a way 
that the data is shared. The BG test results are used to calculate a corrective insulin dose 
automatically, and the insulin delivery device is automatically set for use. Also, carbs 
can be entered manually so that Meal Insulin Boluses are similarly calculated and pre-set. 

40 Digital memory is available in either the meter or the insulin delivery device so that the 
combined BG and insulin history can be downloaded by the Practitioner. 

There are several types of insulin available. The input form contains input boxes for the 
brand names or generic names of the two types of insulin prescribed. 
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The formulas deal with basal amounts in each time interval, rather than basal rates. The 
long time-response of the long-acting insulin makes Basal Rate schedule adjustments 
impractical. 

5 

6.2.1.2.1 BASAL FLOAT 1 

The partnership of Basal Float 1 and Meal Insulin Float 1 may be employed in a similar 
manner to pumps Types D and E. Similarities and differences: As in pumps the 
Practitioner inputs deRxInsl and deMeallns. The input boxes for BR(i) are replaced by a 
10 single box for Basal, a single BR is calculated (one twenty-fourth of Basal) . The values 
of BRrec(i) are shown individually as in pumps. As in pumps, a value is calculated for 
BRaveRec using equation (141). All the BRrx(i) are replaced by a single box for BRrx, 
and Basalrx is calculated (24 times BRrx). After reviewing the recommendations for 
Basal and Basal Rate, the Practitioner inputs BRrx 

15 

6.2.1.2.2 MEAL INSULIN FLOAT 1 

The Meal Insulin Float 1 is similar to that used with Type D pumps. The Enhanced 
insulin nomenclature system is used. 

20 The value of BasalRx from the previous section (or some other source) is used as input 
by the Meal Insulin Float 1 algorithm. 

dBasal(i) = (Basalrx - Basal) * dt(i) / 24 (170) 

25 The present invention then calculates a schedule of recommended CIRrec(i) values, using 
a Meal Insulin Float. The formula is adapted from equation (103) by applying equation 
(132): 

CIRrec(i) = CarbSh(i) / [ eMeallns(i) + (deRxInsl * TBCorIns(i+l)/TBCorTot- 
30 (BaslRx - Basal)*dt(i)/24 ] (171) 

6.2.1.2.3 MEAL INSULIN FLOAT 2 (uses AMCorlns(i) as an error ) 

This algorithm makes use of after-meal BG testing. These tests and associated boluses 
35 divide the interval into two parts. 

Define: Timeb(i)= The time of the after-meal bolus after the ith time boundary. 

dtb(i) = Time(i+1) - Timeb(i) 
dta(i) = Timeb(i) - Time(i) 

40 

The algorithm uses AMCorlns(i) as an "error" indicator for the first part of the interval, 
where the float is calculated. A figure for BasalRx is needed. The Practitioner places it in 
an input box after reviewing the invention's calculated recommendations, which may be 
one or more of the following: 



(172) 
(173) 
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BasalRec may be obtained from the carb-free second part of the interval, which uses 
TBCorIns(i+l) as an error term for dBasal(i) The "Enhanced" insulin terminology 
is not used. In the second part of the interval: 
5 BRrec(i) = BR(i) + Krxlnsl * TBCorIns(i+l) /dtb(i) (174) 

or 

This result from a "reliable" interval may be converted to a full day's basal as 
follows: BasalRec = BRreliable*Kf * 24 (176) 

or 

10 BasalRec may be obtained from the sum of the results above: 

BasalRec = SUM over i ( BRrec(i) * dt(i)) (176) 

or 

BasalRec may be obtained from a Basal Float 1 calculation in a "reliable" interval. 
This result is converted in the same manner as equation (171) 
15 or 

Another s ource o f BasalRec is the total recommended Basal from a Basal Float 1 
algorithm, converted in the same manner as equation (172) 

Whatever the source for BasalRec, the Practitioner inputs BasalRx. The following 
20 statement can be said of the first part of the interval: 

MeallnsRec(i) = Meallns(i) + Krxlnsl *AMCorIns(i) - Basal * dta(i)/24 (177) 

CERrec(i) = CarbSh(i) / MeallnsRec(i) (178) 

The factor Krxlnsl may be obtained from dRxInsl input by the practitioner as described 
25 earlier: 

Krxlnsl = dRxInsl / (AMCorTot +TBCorTot) (1 79) 

Alternatively, if the method of equation (170) is used, Krxlnsl may be replaced by two 
"K-factors", one for each part of the time interval: KrxInsAM and KrxInsTB: 
30 First, Basal Rate is determined from the second part of the interval: 

BRrec(i) = BR(i) + KrxInsTB * TBCorIns(i+l) /dtb(i) (1 80) 

Then Meal Insulin is determined from the first part of the interval: 
35 MeallnsRec(i) = Meallns(i) + KrxInsAM *AMCorIns(i) 

- dta(i) * KrxInsTB * TBCorIns(i+l) /dtb(i) (181) 

CIRrec(i) = CarbSh(i) / MealInsRec(d,i) (1 82) 

40 This alternative method may shift the Basal/TDD ratio and may be useful if such a result 
is intended. 
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6.2.1.3 LIMITED DOMAIN OF deRxInsl and the SAFETY-NET FORMULA for 
MULTIPLE DAYS' DATA VERSIONS 

6.2.1.3.1 LIMITED DOMAIN OF deRxInsl 

The Float formulas using Multiple Days' Data for whole intervals work only within 
5 certain limits in which the changes called for are less than TBCorTot and in the same 
direction (same sign). 

Both these conditions are contained in the requirement that 

TBCorTot / deRxInsl > 1 (1 83) 

10 If the parameters go beyond this limit then nonsensical results occur. To prevent this 
from happening and to provide a further margin of safety, there are input limits placed 
upon the parameters to prevent the Practitioner from inputting numbers out of the 
acceptable range. The input limit is in the form of the parameter KrxInslMax, The limit 
is in the form: 

15 

TBCorTot / deRxInsl >- 1 / KrxInslMax (1 84) 

If KrxInslMax <= 1, then the nonsensical results are avoided. The value of 0.5 is 
currently in use. 

20 

6.2. 1 .3.2 SAFETY-NET FORMULA: 
6.2.1.3.2.1 For Basal Float 1 

25 As mentioned above, the Float formulas using Multiple Days' Data work only within 
certain limits. The first line of defense is the input limits mentioned above. However, 
there MAY BE ways to circumvent the input limits, on purpose or accidentally. To 
protect against this possibility, there is a second line of defense, the safety net formulas. 
Each of the Float formulas actually uses two formulas for deRxInsl(i), each over a 

30 different domain of deRxInsl. 

Small Domain: This is the "acceptable" domain, using the formulas discussed so far. If 
TBCorTot/deRxInsl>l, then the formula below is used. It produces accurate results over 
its limited domain in which the changes called for a re less than or equal to the Total 
35 Enhanced Corrective Insulin and in the same direction. It apportions to each time 
interval a small amount of the deRxInsl proportional to the time interval's share of the 
TBCorTot. Most patient/practitioner interactions are in the small domain. The "small 
domain" formula has been discussed earlier. It is re-written below. The suffix "sd" is 
added to the name the name to show what it is: 

40 

deRxInslsd(i) = deRxInsl*TBCorIns(i+l)/TBCorTot (1 85) 
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Large Domain: The formula introduced below is a backup to the input limits placed upon 
deRxInsl. If the desired change to deRxInsl is not in the small domain, i.e. if 
TBCorTot/deRxInsKl, then the formula below is used. It produces a less accurate but 
safe result when the changes called for are larger than TBCorTot or in the direction 

5 opposite to TBCorTot. The right-hand side starts with the full corrective insulin for the 
time interval Then the "overshoot" of deRxInsl over TBCorTot is apportioned among 
the time intervals, but the apportioning fraction includes Basal as well as Enhanced 
Corrective Insulin, making it a less sensitive fraction, but less exact. In its role as a 
backup to the input limits, the large domain formulas are designed never to be used. The 

10 Safety-Net formula for Type D pumps, Basal Float 1 Algorithm 6.2. 1 . 1 . 1 . 1 , is as follows: 

deRxInslld(i) = TBCorIns(i+l) + (deRxInsl - TBCorTot) *(TBCorIns(i+l) + 
Basal(i))/(TBCorTot + Basal) (186) 

15 The transition between these two formulas is accomplished by a logical variable called 
Swtch, which has a value of zero (0) in the small domain and one (1) in the large domain. 

IF (TBCorTot /deRxInsl>=l) THEN Swtch=0 ELSE Swtch=l (187) 

20 The invention accomplishes this as follows: 

Swtch = IF(TBCorTot /deRxInsl >=1,0,1) (188) 

Swtch helps to put the two equations together into a single equation, which is switched 
25 from (236) to (238) as needed. This is shown below: 

deRxInsl(i)= swtch * deRxInslld(i) + (l-swtch)*deRxInslsd(i) (189) 

After substituting the large and small domain formulas into this, it reduces to the 
30 following equation, which requires fewer computer fields: 

deRxInsl(i) = swtch*TBCorIns(i+l) + (deRxInsl - Swtch*TBCorTot)*(TBCorIns(i+l) + 
Swtch*Basal(i)) / (TBCorTot+Swtch*Basal) (190) 

35 This equation in turn is substituted into one of the Float equations, e.g. equation (140) to 
yield the result below, which produces a schedule of recommended basal rates. 

BRrec(i) = [ Basal(i) + Swtch*TBCorIns(i+l) + (deRxInsl- 
Swtch*TBCorTot)*(TBCorIns(i+l)+Swtch*Basal(i))/(TBCorTot+Swtch*Basal)- 
40 dMealIns*CarbSh(i)/CarbShTot ] /dt(i) (191) 

For Pump Types A, TBCorlnsA(i) is substituted for TBCorlns(i) and TBCorTotA for 
TBCorTot in the above equation. 
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6.2.1.3.2.2 For Meal Insulin Float 1 

The same definitions of the two domains for deRxInsl are used in the Meal Insulin Float 
formulas, Algorithm 6.2. 1 . 1 . 1 .3 Leading to: 

5 

eMeallnsRec(i) = eMeallns(i) + Swtch*TBCorIns(i+l) + ( deRxInsl - Swtch*TBCorTot ) 
*[ TBCorIns(i+l)+ Swtch*eMealIns(i) ] / [ TBCorTot + Swtch*eMealInsTot ] 
- ( BRrx(i) - BR(i) )*dt(i) (192) 

10 This is converted to CIRrec(i) by an adaptation of equation (22): 

CIRrec(i)= CarbSh(i) / eMeallnsRec(i) (193) 

CIRrec(i) = CarbSh(i) / [ eMeallns(i) + swtch* TBCorTot + (deRxInsl - 
15 swtch*TBCorTot)*(TBCorIns(i+l)+ 

swtch*eMealIns(i))/(TBCorTot+swtch*eMealInsTot) - dt(i)*(BRrx(i) 

-BR(i))] (194) 

Where Swtch is the same "switching parameter" introduced earlier by equation (183) to 
20 shift the formula between two different domains. 



6.2.1.3.2.3 For Meal Insulin Float 2 (uses AMCorlns(i) as error ) 
Small Domain: 
25 Over the whole interval: 

dRxInslsd(i) = [dRxInsl / (AMCorTot +TBCorTot)] * (AMCorlns(i) + 

TBCorlns(i)) (195) 
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This can be thought of as an equation in two parts: 
For the last part of the interval: 

dRxInslsdb(i) = [dRxInsl / (AMCorTot +TBCorTot)] * TBCorlns(i) = 

dBasalb(i) (196) 

35 Note that this is all Basal Insulin. 

For the first part of the interval: 

dRxInslsda(i) = [dRxInsl / (AMCorTot +TBCorTot)] * AMCorlns(i) ( 1 97) 
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Large Domain: 

For the last part of the interval: 

dRxInslldb(i) = TBCorlns(i) +(dRxInsl - AMCorTot -TBCorTot)* 
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( TBCorIns(i)+ BR(i)*dtb(i))/(AMCorTot +TBCorTot + Basal) 



(198) 



For the first part of the interval: 

dRxInsllda(i) = AMCorlns(i) +(dRxInsl - AMCorTot -TBCorTot)* 
5 ( AMCorlns(i) + BR(i)*dta(i))/(AMCorTot +TBCorTot + Basal) (199) 

The sum of the two corrective insulin totals must be used in the switch parameter. 
Swtchm = IF((AMCorTot + TBCorTot) /dRxInsl >=1 ,0, 1 ) (200) 

10 

6.2.2 AUTOMATIC MULTIPLE DAYS' DATA (DIGITAL ADVISOR) 
The descriptions so far in Version 2 describe the "manual" version of the Multiple Past 
Days' algorithm. The automatic features described below are to remove the pauses for 
1 5 human judgement and input. 

6.2.2.1 AUTOMATION OF deRxInsl: 

To automate the parameter deRxInsl, the present invention makes use of the fact that it is 
more difficult to control the diabetes of a patient whose BG tests, have a high percent 
20 standard deviation, (BGsd/BGmean) compared to the database norm. The present 
invention uses a simple ramp function for this as shown below: 

Some new parameters: 

Fins Auto is a multiplier for use in the formula for deRxInsAuto as follows: 
25 deRxInsAuto=FInsAuto * KrxInslMax*TBCorTot (201) 

where 

KrxInslMax is a positive constant <=1, introduced earlier. It is subject to change by 
30 the programmers if necessary. The present value is 0.5. 

Two definitions: 

PmPctBGsd: Database "population" mean of the quantity (BGsd/BGmean) . 
35 PsdPctBGsd: Database "population" standard deviation of the quantity (BGsd/BGmean), 
i.e. the database standard deviation of the personal standard deviations. 

In the present invention, FInsAuto is a ramp function as follows: 

40 IF BGsd/BGmean < PmPctBGsd + PsdPctBGsd then FInsAuto = 1 

IF PmPctBGsd + PsdPctBGsd <= BGsd/Bgmean < =PmPctBGsd + 2 PsdPctBGsd, 

THEN FInsAuto = [ 1 - (BGsd/Bgmean - PmPctBGsd - PsdPctBGsd) / PsdPctBGsd ] 
IF BGsd/BGmean > PmPctBGsd + 2 PsdPctBGsd then FInsAuto = 0 (202) 
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This function describes a flat region followed by a ramp down to zero 
as shown in Figure 3. 



6.2.2.2 AUTOMATION OF deMeallns: 
5 The present invention makes use of the parameter dBslToTgt, which is a general advisory 
parameter that shows what change in Basal is necessary to achieve BoTTgt. 

dBslToTgt - TDD * BoTTgt - Basal (203) 

10 Two steps are employed to automate dMeallns: 

The first equation below assigns to the change an absolute value equal to the minimum 
absolute value of deRxInsl or dBaslToTgt. The reason is to avoid overshooting either 
BoTTgt or deRxInsl. The second equation uses the result from the first, together with 
15 equation (57), to calculate a max value for deMeallns Auto; then it limits deMeallns Auto 
to Kmauto times its max value if it is positive (in the direction of adding insulin). The 
parameter Kmauto is subject to adjustment by the programmers. 

dBasalAuto = IF[ABS(dBaslToTgt/deRxInsl) <1, ABS(dBaslToTgt), 
20 ABS(deRxInsl)]*sgn(dBaslToTgt) (204) 

deMealInsAuto=IF[(deRxInsl>dBasalAuto,Kmauto*(deRxInsl-dBasalAuto), (deRxInsl- 
dBasalAuto)] ' (205) 

25 where Kmauto is a positive constant < =1 , presently set at 1 . 

The value of BoTTgt may be obtained from the AIM study's value for Kb or from 
BoTTgtRec from Section 6. 1 .2. 1 Basal-to-Total Ratio 

30 6.2.2.3 AUTOMATIC ROUNDING OF CIR FOR PUMP TYPE A: 

In the manual mode for Type A, a two-step process is required to obtain a half-integer 
value for CIR in the range below 8. The automatic feature accomplishes this by 
multiplying the floating point number dCIRcalcA by two, then rounding off to a whole 
number, and then dividing by two. The result is input into the reverse equation to correct 

35 the dMeallns for the change. 

6.3 NON-FLOAT ALGORITHM 

This algorithm uses the BoTFbk(d) factor to determine both Basal Rate and CIR. This 
40 type of algorithm is shown below in a Daily Update Version: 

eMealIns(d,i) = eMeallns(d-U) + KrxInsl*(l-BoTFbk(d-l))* 

TBCorIns(d-l,i+l) (206) 
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CIR(d,i) = CIR(d-l,i) - KrxInsl*(l-BoTFbk(d-l)) * eCorIns(d-l,i+l) * CIR(d- 
l,i)/eMealIns(d-l,i) (207) 



BR(d,i) = BR(d-l,i) + KrxInsl*BoTFbk(d-l)*TBCorIns(d-l,i+l) (208) 

5 

This method has the advantage of being very accurate, but has the disadvantage of 
assigning large changes in meal insulin to time intervals that have little to no carb 
consumption. For this reason, the version is suitable only for the intervals containing the 
main meals. 

10 

7. TIME INTERVALS: 

Theoretically, there should be no limit on the number of time intervals in the day. 
However, for reasons of screen space, there is a limit on the number of time intervals that 
may be used. Also, for an individual patient, fewer than the screen's maximum may be 
15 used. The invention needs to find the last one. The last time boundary of the day (before 
midnight) is called Tmax. It is found by a routine that picks out the time boundary with 
the highest 24-hour clock time. This routine is incorporated into a single formula of 
nested IF statements shown below, (note: In MS Access the NZ function makes nulls 
and zeros act the same): 

20 

Tmax= 

nf(NZ([Timel])>NZ([Time2]),[tt^ 
Time3])>NZ([Time4]),[Time3^ 

>NZ([Time6]),[Time5],nf(NZ([Time6])>NZ([Time7]),[Time6],[Time7])))))) (209) 

25 

The Time Intervals, dt(i), are found by subtracting the consecutive time boundaries. 

dt(i) = time(i+ 1 ) - time(i) copy of (4) 

Special mention should be made of the interval surrounding midnight. 
30 dtO -24 + timel - Tmax copy of (5) 



7.1 FINDING THE DAY'S LAST ENTRY 

35 Since people rarely go to bed exactly at midnight, it often occurs that some parameters 
normally associated with bedtime affect the early morning time interval. For instance, 
Meal Insulin taken at bedtime affects the Corrective Insulin of the mid-sleep time 
boundary. The formula below will find the last entry of the day: 

40 [MeallnsLast] = IIF([Time2]=[Tmax], [MealIns2],0)+IIF([Time3]=[Tmax],[MealIns3],0) 
+IIF([Time4]=[Tmax],[MealIns4],0) + IIF([Time5]=[Tmax],[MealIns5],0) + 
IIF([Time6]=[Tmax],[MealIns6],0) + IIF([Time7]=[Tmax],[MealIns7],0) (210) 
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The recommended values for the last time interval in the day need to be re-located on the 
screen. The following is an example using Recommended Basal Rate, BR: 

BRrec(i) = HFflTime(i)] = [Tmax], [BRrecLast], [BRrec(i) regular formula]) (211) 

5 

Other alternative embodiments will become apparent to those skilled in the art to which 
an exemplary embodiment pertains without departing from its scope and spirit. 
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ADDENDUM A 



SQL listing for Multiple Day's Data Version for Insulin Pumps (Type C). This Version 
of the Invention is embodied in a Microsoft Access database for use by the Practitioner as 
5 he/she interacts with a patient. 



SELECT TI.Tp_PatientID, TP.SSN, TP.Name, TP.Endo, TP.DOB, TP.Sex, 
TP.OnsetDMdate, TP.PumpStartDate, TP.PhoneNos, TP.Email, 

10 

TI.IntractionID, TI.IntrDate, 

TI.Hgt, Tl.Wt, TI.HbAlc, Tl.Creatinine, TI.Microalb, 
TI.BGmean, TI.BGsd, TI.BGpd, TI.NBGs, 

1 5 TI.TDDa, TI.TDDb, TI.TDDc, TI.TDDd, TI.TDDe, TI.TDDf, 

(NZ([TDDa])+NZ([TDDb])+NZ([TDDc])+NZ([TDDd])+NZ([TDDe])+NZ([TDDf]))/(- 
IsNumeric([TDDa])-IsNumeric([TDDb])-IsNumeric([TDDc])-IsNumeric([TDDd])- 
IsNumeric([TDDe])-IsNumeric([TDDfJ)) AS TDDavg, 

20 [Basal]/[TDDavg] AS BasalOvrTDD, 
0.48* [TDDavg] AS BasaLADvI, 
[BasalAEvi]-[Basal] AS dBslToAIM, 
TI.Basalrx, 

TI.TargetBG, TI.TargetBGrx, 
25 TI.CF, 1708/[TDDavg] AS CFaim, TI.CFrx, 
TI.CIR, 2.81*[wt]/[TDDavg] AS CIRaim, 

IIf([sex]="m",(([Hgt]-60)*6+l 06)* 1 3/8,(([Hgf]-60)*5+l 00)* 1 3/8) AS [Carbspd(Hgt)] , 
l-[CarbShTot]*0.8/[Carbspd(Hgt)] AS BoTTgtRec, 
[CarbShTot]/[eMealInsTot] AS [CIR(carb)], TI.MeallnsTotAPoc, 
30 [CarbShTot]/[MealInsTotARx] AS [CIR(MealInslRx)], 
TLTypelnsShort, TI.TypelnsLong, 
TI.Diet, TI.PumpType, Tl.iComment, 

TI.TimeLabl2, TI.TimeLabl3, TI.TimeLabl4, TI.TimeLabl5, TI.TimeLabl6, 
35 TI.TimeLabl7, TI.TimeLabl8, 

0 AS Timel, TI.Time2, TI.Time3, TI.Time4, TI.Time5, TI.Time6, TI.Time7, TI.Time8, 

IIf(NZ([Time2])>NZ([Time3]) ) [Time2],nf(NZ([Time3])>NZ([Time4]),[Time3],IIf(NZ([ 
40 Time4])>NZ([Time5]),[Time4],IIf(NZ([Time5])>NZ([Time6]),[Time5] ) IIf(NZ([Time6]) 
>NZ([Time7]),[Time6],nf(NZ([Time7])>NZ([Time8]) > [Time7],[Time8])))))) ASTmax, 

TI.CarbSh2, TI.CarbSh3, TI.CarbSh4, TI.CarbSh5, TI.CarbSh6, TI.CarbSh7, 
TI.CarbSh8, 
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NZ([CarbSh2])+NZ([CarbSh3])+NZ([CarbSh4])+NZ([CarbSh5])+NZ([CarbSh6])+NZ([ 
CarbSh7])+NZ([CarbSh8]) AS CarbShTot, 

IIf([Time2]=[Tmax],[CarbSh2],0)+nf([Time3]=[Tmax] ) [CarbSh3],0)+IIf([Time4]=[Tma 
x],[Carbsh4],0)+IIf([Time5]=[Tmax],[CarbSh5],0)+nf([Time6]=[Tmax],[CarbSh6],0)+II 
f([Time7]=[Tmax],[CarbSh7],0)+NZ([CarbSh8]) AS CarbShLast, 

TI.MealIns2, TI.MealIns3, TI.MealIns4, TI.MealIns5, TI.MealIns6, TI.MealIns7, 
TI.MealIns8, 

TI.AMCorIns2, TI.AMCorIns3, TI.AMCorIns4, TI.AMCorIns5, TI.AMCorIns6, 

TI.AMCorIns7, 

TI.AMCorIns8, 

TI.TBCorIns2, TI.TBCorIns3, TI.TBCorIns4, TI.TBCorIns5, TI.TBCorIns6, 
TI.TBCorIns7, TI.TBCorIns8, 

NZ([TBCorIns2])+NZ([TBCorIns3])+NZ([TBCorIns4])+NZ([TBCorIns5])+NZ([TBCor 
Ins6])+NZ([TBCorIns7])+NZ([TBCorIns8]) AS TBCorTot, 



[MealIns2]+[AMCorIns2] AS eMealIns2 
[MealIns3]+[AMCorIns3] AS eMealIns3 
[MealIns4]+[AMCorIns4] AS eMealIns4 
[MealIns5]+[AMCorIns5] AS eMealIns5 
[MealIns6]+[AMCorIns6] AS eMealIns6 
[MealIns7]+[AMCorIns7] AS eMealIns7 
[MealIns8]+[AMCorIns8] AS eMealIns8 

NZ([eMealIns2])+NZ([eMealIns3])+NZ([eMealIns4])+NZ([eMealIns5])+NZ([eMealIns 
6])+NZ([eMealIns7])+NZ([eMealIns8]) AS eMeallnsTot, 

IIf([Time2]=[Tmax],NZ([eMealIns2]),0)+IIf([Time3]=[Tmax],NZ([eMealIns3]),0)+IIf([ 
Time4]=[Tmax],NZ([eMealIns4]),0)+IIf([Time5]=[Tmax],NZ([eMealIns5]),0)+IIf([Tim 
e6]=[Tmax],NZ([eMealIns6]),0)+IIf([Time7]=[Tmax],NZ([eMealIns7]),0)+NZ([eMealIn 
s8]) AS eMeallnsLast, 

[Time2]+24-[tmax] AS dtO, 
IIf(NZ([Time2])>0,[time2]-[timel],0) AS dtl, 
IIf(NZ([Time3])>0 > [time3]-[time2],0) AS dt2, 
IIf(NZ([Time4])>0 > [time4]-[time3],0) AS dt3, 
IIf(NZ([Time5])>0,[time5]-[time4],0) AS dt4, 
IlfCNZCtTimedJ^O.ttimeei-ttimeS]^) AS dt5, 
IIf(NZ([Time7])>0,[time7]-[time6],0) AS dt6, 
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IIf(NZ([Time8])>0,[time8]-[time7],0) AS dt7, 

TI.BR1, TI.BR2, TI.BR3, TI.BR4, TI.BR5, TI.BR6, TI.BR7, TI.BR8, 

[BRl]*[dtO] AS BaslO, 
[BR2]*[dt2] AS Basl2, 
[BR3]*[dt3] AS Basl3, 
[BR4]*[dt4] AS Basl4, 
[BR5]*[dt5] AS Basl5, 
[BR6]*[dt6] AS Basl6, 
[BR7]*[dt7] AS Basl7, 

NZ([Basl7])+NZ([Basl6])+NZ([Basl5])+NZ([Basl4])+NZ([Basl3])+NZ([Basl2])+NZ([B 
aslO]) AS Basal, 

TI.BG2, TI.BG3, TI.BG4, TI.BG5, TI.BG6, TI.BG7, TI.BG8, 

TI.CIR1, TI.CIR2, TI.CIR3, TI.CIR4, TI.CIR5, TI.CIR6, TI.CIR7, TLCIR8, 

IIf([Time2]=[Tmax],[CIR2],0)+IIf([Time3]=[Tmax],[CIR3],0)+IIf([Time4]=[Tmax],[CI 
R4],0)+IIf([Time5]=[Tmax],[ClR5] ) 0)+IIf([Time6]=[Tmax],[CIR6],0)+IIf([Time7]=[Tm 
ax],[CIR7],0)+NZ([CIR8]) AS CIRLast, 



TI.deRxInsl, 
TI.deMeallns, 

[deRxInsl]-[deMealIns] AS dBasal, 
[CarbShTot]/([eMealInsTot]+[deMealIns]) AS CIRrecC, 
TI.CIRrx, 

IIf(0<=[deRxInsl]/[CorTot]<=l,0,l) AS swtch, 

([BaslO]+[swtch]*[TBCorIns2]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns2]+[swtch]*[BaslO]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*[eMealInsLast]/[eMealInsTot])/[dtO] AS 
BRrecAO, 

[BRrecAO] AS BRrecAl, 

([Basl2]+[swtch]*[TBCorIns3]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns3]+[swtch]*[Basl2]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns2])/[eMealInsTot])/[dt2] 
AS BRrecA2, 
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([Basl3]+[swtch]*[TBCorIns4]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns4]+[swtch]*[Basl3]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns3])/[eMealInsTot])/[dt3] 
AS BRrecA3, 

([Basl4]+[swtch]*[TBCorIns5]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns5]+[swtch]*[Basl4]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns4])/[eMealInsTot])/[dt4] 
AS BRrecA4, 

([Basl5]+[swtch]*[TBCorIns6]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns6]+[swtch]*[Basl5]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns5])/[eMealInsTot])/[dt5] 
AS BRrecA5, 

([Basl6]+[swtch]*[TBCorIns7]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns7]+[swtch]*[Basl6]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns6])/[eMealInsTot])/[dt6] 
AS BRrecA6, 

([Basl7]+[swtch]*[TBCorIns8]+([deRxInsl]- 
[swtch]*[TBCorTot])*([TBCorIns8]+[swtch]*[Basl7]) 

/([TBCorTot]+[swtch]*[BaslTot])-[deMealIns]*NZ([eMealIns7])/[eMealInsTot])/[dt7] 
AS BRrecA7, 

TI.BRrxl, TI.BRrx2, TI.BRrx3, TI.BRrx4, TI.BRrx5, TI.BRrx6, TI.BRrx7, TI.BRrx8, 
0 ASTimeRxl, 

TI.TimeRx2, TI.TimeRx3, TI.TimeRx4, TI.TimeRx5, TI.TimeRx6, TI.TimeRx7, 
TI.TimeRx8, 

IIf(NZ([TimeRx2])>NZ([TimeRx3]),[TimeRx2],nf(NZ([TimeRx3])>NZ([TimeRx4]),[Ti 
meRx3]JIf(NZ([TimeRx4])>NZ([TimeRx5]),[TimeRx4],IIf(NZ([TimeRx5])>NZ([Time 
Rx6]),[TimeRx5],nf(NZ([TimeRx6])>NZ([TimeRx7]),[TimeRx6],IIf(NZ([TimeRx7])> 
NZ([TimeRx8]),[TimeRx7],[TimeRx8])))))) AS TRxmax, 

IIf(NZ([TimeRx2])>0,[TimeRx2]+24-[Trxmax],0) AS dTrxO, 
IIf(NZ([TimeRx2])>0,[TimeRx2]-[TimeRxl],0) AS dTrxl, 
IIf(NZ([TimeRx3])>0,[TimeRx3]-[TimeRx2],0) AS dTrx2, 
IIf(NZ([TimeRx4])>0,[TimeRx4]-[TimeRx3],0) AS dTrx3, 
IIf(NZ([TimeRx5])>0,[TimeRx5]-[TimeRx4],0) AS dTrx4, 
IIf(NZ([TimeRx6])>0,[TimeRx6]-[TimeRx5],0) AS dTrx5, 
IIf(NZ([TimeRx7])>0,[TimeRx7]-[TimeRx6],0) AS dTrx6, 
IIf(NZ([TimeRx8])>0,[TimeRx8]-[TimeRx7],0) AS dTrx7, 
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[CarbShLast]/([eMealInsLast]+[deRxInsl]*[TBCorIns2]/[TBCorTot]-([BRrxl]- 
[BRl])*[dtO]) AS CIRrecO, 

[CarbSh2]/([eMealIns2]+[deRxInsl]*[TBCorIns3]/[TBCorTot]-([BRrx2]-[BR2])*[dt2]) 
AS CIRrec2, [CarbSh3]/([eMealIns3]+[deRxInsl]*[TBCorIns4]/[TBCorTot]-([BRrx3]- 
[BR3])*[dt3]) AS CIRrec3, 

[CarbSh4]/([eMealIns4]+[deRxInsl]*[TBCorIns5]/[TBCorTot]-([BRrx4]-[BR4])*[dt4]) 
AS CIRrec4, [CarbSh5]/([eMealIns5]+[deRxInsl]*[TBCorIns6]/[TBCorTot]-([BRrx5]- 
[BR5])*[dt5]) AS CIRrec5, 

[CarbSh6]/([eMealIns6]+[deRxInsl]*[TBCorIns7]/[TBCorTot]-([BRrx6]-[BR6])*[dt6]) 
AS CIRrec6, [CarbSh7]/([eMealIns7]+[deRxInsl]*[TBCorIns8]/[TBCorTot]-([BRrx7]- 
[BR7])*[dt7]) AS CIRrec7, 

TI.eMeallnsRxl, TI.eMealInsRx2, TI.eMealInsRx3, TI.eMealInsRx4, TI.eMealInsRx5, 
TI.eMealInsRx6, TI.eMealInsRx7, 

TLCIRixl, TI.CIRrx2, TI.CIRrx3, TI.CIRrx4, TI.CIRrx5, TI.CIRrx6, TI.CIRrx7, 
TI.CIRrx8, TI.InptFormType 

FROM TP INNER JOIN TI ON TP.PatientID = TI.Tp_PatientID 
ORDER BY TI.IntrDate DESC; 
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